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ABSTRACT 


The  Camegie-Mellon  University  team  has  completed  the  second  two 
year  program  in  advanced  composites  technology.  The  program  has  had 
significant  impact  in  several  areas  as  the  CMU  team,  working  closely  with 
the  engineers  at  General  Dynamics,  Convair  Aerospace  Division,  Fort 
Worth,  has  continued  work  on  several  projects.  The  work  reported  herein 
included:  continuing  development  of  university-industry  interaction; 
development  of  educational  material  including  reports  on  manufacturing 
methods,  Weibull  statistics,  and  stresses  due  to  an  elliptical  hole  in 
an  anisotropic  plate;  and,  advanced  research  into  synthesis  procedures 
for  mechanically  fastened  joints  and  elastic  fracture  mechanics  for 
laminated  composite  plates.  The  research  projects  demonstaate  useful 
design  capabilities,  new  analystic  and  numerical  results,  and  experiment¬ 
al  data. 
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CHAPTER  I 


REVIEW  OF  THE  INTERACTIVE  PROGRAM 

1.1  Introduction 

This  final  report  is  divided  into  three  parts:  Volume  I  reviews 
some  aspects  of  the  Interactive  Program  during  the  contract  period; 
Volumes  II  and  III  are  Ph.D.  theses  written  during  the  course  of  the 

"k 

current  and  preceding  contract  period  [1].  The  following  chapters 
in  this  volume  consist  of  selected  student  reports  from  various 
project  activities. 

The  Interactive  Program  has  been  a  successful  program  of  coupling 
faculty  and  students  at  Carnegi e-Mel  Ion  University  (CMU),  The  Air 
Force  Materials  Laboratory,  and  industry,  most  notably  General  Dynamics 
Corporation,  Convair  Aerospace  Division,  Fort  Worth.  The  program  has 
focused  on  the  application  of  mechanics  capabilities  at  CMU  to  the 
stress  and  strength  analysis  of  advanced  composite  structures.  The 
broad  objectives  of  the  program  have  been  the  following: 

1.  Creation  of  effective  means  of  communication  between  CMU  and 
industry; 

2.  Involvement  of  the  CMU  team  of  students  and  faculty  in  the 
definition  and  solution  of  fundamental  problems  arising  from 
the  application  of  advanced  composites  in  aerospace  structures; 
and 


*References  are  denoted  by  brackets  and  are  located  at  the  end 
of  each  volume. 
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3.  Development  by  the  CMU  team  of  new  analysis  capabilities  and 
results,  strength  criteria,  design  mdthods,  and  educational 
material  for  advanced  composites  technology. 

Two  elements  were  keys  to  the  success  of  the  CMU  program.  The 
first  was  the  development  of  an  educational  basis  for  the  student 
project  activities.  Due  to  the  relatively  new  nature  of  composites 
technology,  this  was  accomplished  as  an  early  part  of  the  preceeding 
contract.  The  basic  outline  of  the  educational  material  and  the  de¬ 
scription  of  the  program  organization  and  timing  are  reviewed  in  [1]. 

The  second  element  for  success  was  the  development  of  a  record 
of  success  and  expertise,  coupled  with  a  level  of  external  visibility 
sufficient  to  give  credibility  to  the  CMU  team.  The  major  contribu¬ 
tion  for  this  success  has  been  the  use  of  a  major  portion  (20-25%)  of 
the  program  budget  for  travel.  These  trips  generally  involved  the 
students  and/or  the  principal  investigator  travelling  to  industry. 

By  making  frequent  visits  to  industry  and  program  review  meetings  it 
was  possible  to  develop  contacts,  credibility  and  visibility  necessary 
to  the  program's  success.  Further  aspects  of  the  program  visibility 
are  described  below. 

The  CMU-developed  philosophy  also  contributed  materially  to 
the  success  of  the  program.  It  was  initially  decided  to  concentrate  on 
problems  which  had  a  clear  engineering  significance  as  opposed  to 
more  "academic"  pursuits.  Thus,  initially  at  least  post  graduate  students 
were  not  wedded  to  a  particular  research  program.  Further,  it  was 
decided  to  maintain,  to  the  extent  possible,  a  framework  of  project 
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activities  closely  attuned  to  current  CMU  faculty  capabilities. 

These  included  fracture  mechanics  and  elastic  stress  analysis.  Ample 
project  problems  were  found  that  were  based  in  these  areas  but  which 
had  immediate  pay-off  potential  for  industry. 

The  two  most  significant  derivatives  of  this  philosophy  are  attached 
as  separate  volumes  of  this  report,  both  of  which  are  Ph.D.  theses 
written  as  a  result  of  projects  started  four  years  ago.  The  first 
develops  an  efficient  synthesis  procedure  for  mechanical  joints  in 
composites;  the  second  concerns  a  macromechanics  basis  for  an  engineering 
fracture  mechanics  of  advanced  composites. 

1.2  Research  Results  of  the  Program 

As  noted  above,  the  CMU  team  undertook  a  number  of  projects  in 
the  general  area  of  strength  analysis  of  advanced  composites.  The 
approach  used  was  to  retain  the  simple  basis  of  two  dimensional 
elasticity  and  elastic  fracture  mechanics.  Several  papers  [2-6] 
were  published  during  the  contract  to  provide  visibility  for  the 
CMU  efforts.  While  much  of  this  work  retains  a  speculative  flavor, 
particularly  on  the  question  of  fracture  criterion,  it  has  received 
some  positive  response  from  industry.  As  an  example.  General 
Dynamics  engineers  have  made  use  of  fracture  mechanics  methodology  to 
define  a  fail-safe  composite  panel. 

The  research  results  have  been  described  in  several  technical 
formats.  During  the  current  report  period  papers  [2,6,7]  were  pre¬ 
sented  at  AIAA/ASME  Structures,  Structural  Dynamics,  and  Materials 
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conferences  in  1971,  1972,  and  1973.  One  paper  [3]  was  presented  at 
the  1972  ASTM  meeting  on  Analysis  of  Test  Methods  for  Advanced  Fibers  and 
Their  Composites.  These  research  results  were  also  discussed  at  numerous 
Air  Force  review  meetings  and  technical  committee  meetings. 

In  October  1972  Drs.  Cruse,  Swedlow,  and  Halpin  organized  the 
first  Colloquium  on  Structural  Reliability  [9],  under  partial 
support  of  the  Air  Force  Office  of  Scientific  Research.  This 
meeting  brought  together  experts  in  composites,  fracture  mechanics, 
statistics,  and  structural  reliability  for  three  days  of  valuable 
technical  interaction.  This  meeting  was  held  at  Carnegi e-Mel  Ion 
University.  Currently  about  120  copies  of  the  Pvooeedings  [9]  are 
in  circulation. 

1.3  Conclusions 

It  is  not  entirely  possible  to  define  the  success  of  the  Inter¬ 
active  Program  in  Quantitative  terms.  However  it  was  clearly  shown 
that  it  is  possible  to  create  valuable  interaction  between  the 
university  and  industry  environments,  around  a  well  defined  project 
area. 

Further,  the  program  has  contributed  to  the  technical  knowledge 
in  the  area  of  strength  of  advanced  composite  materials.  This  has 
been  accomplished  through  the  technical  literature  as  well  as  through 
various  technical  meetings.  The  fact  that  certain  people  pick  us 
out  to  disagree  with  is  indicative  of  some  degree  of  success. 
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Finally,  the  program  has  seen  to  members  of  the  CMU  team  placed 
in  the  technical  community  concerned  with  composites.  Dr.  H.  J. 
Konish,  Jr.  is  a  post-doctoral  fellow  at  the  Air  Force  Materials 
Laboratory.  Dr.  J.  P.  Waszczak  has  taken  a  position  with  General 
Dyaamics  in  San  Diego,  California.  The  fact  that  both  individuals 
received  several  offers  attests  to  the  success  of  the  engineering 
emphasis  of  the  Interactive  Program. 
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CHAPTER  II 


FABRICATION  PRACTICES  AND  PROBLEMS  FOR  BORON  AND  GRAPHITE  EPOXY 

2.1  Introduction 

Advanced  composite  materials,  particularly  boron  and  graphite 
epoxy  have  reached  the  stage  in  their  development  where  the  possibili¬ 
ties  and  problems  associated  with  fabrication  of  components  from 
these  materials  are  being  fully  investigated.  Now  that  the  materials 
have  reached  a  satisfactory  level  of  design  reliability,  the  time 
has  come  to  determine  whether  structures  can  be  manufactured  from 
boron  or  graphite  fiber  composites  for  reasonable  costs.  The  key 
to  the  widespread  acceptance  of  these  new  materials  is  going  to 
be  whether  the  additional  performance  gained  by  using  them  is 
worth  the  additional  and  sometimes  unique  problems  of  their  fab¬ 
rication.  This  report  examines  the  major  areas  of  advanced 
composite  fabrication  pointing  out  some  of  the  current  techniques 
employed  by  those  in  the  "industry".  The  major  problems  in  each 
area  are  outlined  as  well  as  advances  by  individual  companies  in 
some  particular  aspect  of  that  area. 
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2.2  Tooling 


The  most  important  parameter  in  the  choice  of  a  tooling  material 
is  the  coefficient  of  thermal  expansion  (a).  The  value  of  a  for 
the  tool  should  be  as  close  as  possible  to  the  coefficient  of  thermal 
expansion  for  the  material  that  is  to  be  cured  on  the  tool  for  a 
number  of  reasons.  First  of  all,  an  improper  match  of  these 
coefficients  could  cause  the  curing  material  to  separate  from  the 
tool  surface  and  therefore  not  cure  with  the  proper  contour.  Secondly, 
the  curing  process  involves  the  bleeding  of  resin  from  the  part  and 
the  tool  often  includes  dams  along  the  perimeter  of  the  piece  to 
prevent  excessive  resin  flow.  If  the  tool  has  a  higher  value  of  a 
than  the  part,  the  dams  may  separate  from  the  perimeter  and  allow 
too  much  resin  to  flow  from  the  matrix.  If  the  value  of  a  is  too 
low  in  the  tool,  the  part  will  try  to  expand  more  than  the  dams 
will  permit,  causing  the  part  to  buckle.  Finally,  a  high  coefficient 
of  thermal  expansion  may  cause  the  tool  to  fail  due  to  fatigue 
caused  by  thermal  stresses  in  high-volume  production. 

The  use  of  conventional  tooling  methods  and  materials  has 
thus  far  proven  to  be  adequate  in  the  manufacture  of  advanced 
composite  structures.  Mild  steel  designated  M-2,  is  one  of  the 
more  popular  materials  since  it  is  relatively  inexpensive  and  has 
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a  value  of  a  reasonably  close  to  that  of  boron  epoxy.  However,  most 
of  the  fabrication  of  composites  thus  far  has  been  done  with  relatively 
smooth,  flat  surfaces.  It  is  possible  that  the  manufacture  of 
components  with  a  high  degree  of  curvature  or  a  number  of  complex 
contours  may  require  tools  with  an  almost  perfect  thermal  match. 

At  least  two  companies  have  undertaken  programs  toward  this 
end.  Boeing/Seattle  (B/S)  is  currently  evaluating  the  use  of  a 
ceramic  castable  compound  which  they  have  named  Comptite  II.  The 
exact  composition  of  Comptite  II  is  Boeing  proprietary  information 
but  it  is  basically  composed  of  a  betaeucryptite  and  a  fused  silica 
binder.  By  varying  the  percentages  of  betaeucryptite,  it  is  possible 
to  vary  the  coefficient  of  thermal  expansion  of  the  material.  In 
this  way,  a  for  Comptite  II  can  be  best  matched  to  a  for  boron  or 
graphite  epoxy.  B/S  reports  that  an  exact  correlation  of  a's  may 
cause  distortion  and  that  the  tool  works  best  when  its  coefficient 
of  thermal  expansion  is  slightly  higher  than  that  of  the  boron 
or  graphite  epoxy  which  is  being  cured.  Comptite  II  is  easy  to 
form  and  demonstrates  high  fatigue  life. 

Another  major  exception  to  steel  tooling  was  performed  by 
General  Dynamics  Convair  Division  in  Fort  Worth  (GD).  In  the 
prototype  manufacture  of  a  graphite  epoxy  F-5  fuselage  component, 

GD  used  a  female  graphite  epoxy  tool  to  exactly  match  coefficients 
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of  thermal  expansion.  Although  tool  costs  were  estimated  at 
two-thirds  the  cost  of  a  comparable  steel  tool,  metal  would  be 
used  in  a  high-volume  production  run  because  it  is  believed  that 
steel  could  withstand  more  cycles  in  the  autoclave. 
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2.3  Lay-Up 

The  crux  of  any  fabrication  technique  is  the  method  by  which 
the  raw  material  is  transformed  into  a  final  shape.  In  the  area 
of  composite  structures,  this  stage  of  operation  is  termed  lay-up 
since  a  component  is  formed  by  the  successive  buildup  of  layers 
of  raw  material. 

2.3.1  VoMcxlption  HateJiiaL 

Boron  and  graphite  epoxy  come  in  two  basic  modes:  3-inch 
tape  on  spools,  by  far  the  more  common  of  the  two,  or  "wide  goods" 
which  vary  in  size  and  shape.  The  boron  or  graphite  fibers  are 
all  oriented  in  a  single  direction  in  the  resin  matrix  and  it  is 
the  varying  angles  of  orientation  of  the  layers  which  gives 
the  material  its  desirable  strength  characteristics  after  curing. 

Boron  and  graphite  epoxy  have  shelf-life  limitations  which 
are  highly  dependent  on  temperature.  Most  manufacturers  specify 
a  3-month  shelf-life  for  tape  stored  at  0°F  or  lower.  Once  the 
boron  or  graphite  reaches  room  temperature,  about  70°F,  it  should 
be  used  within  the  next  10  days.  These  rather  specific  time- 
temperature  requirements  necessitate  the  careful  management  of 
inventories  and  lay-up  procedures  to  ensure  that  the  tape  retains 
its  mechanical  properties.  Some  alternative  methods  for  laying-up 
composite  structures  from  3-inch  tape  are  described  below. 
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2.3.2  GmtMl  Vyncmici  Method 
2.3.2./  Kul-omatdd  McLakim  Lay-Up 

GD  has  practical  experience  in  manufacturing  pieces  with  an 
automated  tape-laying  machine.  The  F-111  horizontal  tail  was 
done  on  a  computer  programmed  basis  using  Narmco  5505  boron 
epo)^.  The  machine  is  preset  to  lay  the  tape  to  the  specified 
tolerances  one  layer  at  a  time.  The  head  is  rotated  for  the  correct 
orientation  of  each  layer  and  the  machine  proceeds  to  lay  tape  on 
tape  until  the  component  is  built  up  to  the  correct  number  of  layers. 

A  typical  tape-laying  machine  head  is  shown  in  Figure  1.  Incorporated 
in  some  models  is  a  light-sensitive  inspection  device  which  detects 
broken  and  crossed  fibers  or  excessive  gaps  between  the  fibers  in 
the  3"  tape.  This  section  of  tape  is  automatically  rejected  and 
the  machine  does  not  lay  any  tape  until  the  section  passing  through 
the  detector  again  meets  the  standards.  This  type  of  machine  is 
highly  sophisticated  and  most  of  the  work  done  by  it  has  been 
purely  experimental  rather  than  actual  production.  GD  believes 
that  the  tape  laying  machine  can  have  a  valuable  application  if  it 
is  used  to  manufacture  many  parts  in  a  single  operation.  For 
example,  GD  has  plans  to  manufacture  the  wing  rib  shown  in  Figure 
2  in  the  following  manner: 


11 


Skin 

1.  Lay-up  a  solid  laminate  30'  x  10'  and  cure  using  the 
vacuum  bag  process. 

2.  Cut  the  large  laminate  into  multiple  panels  of  the 
correct  size. 

(Note):  The  30'  x  10'  laminate  may  be  cut  before  or  after 
curing  depending  on  the  availability  of  composite 
machining  equipment. 

Hat  StiiimoM 

1.  Fabricate  graphite  flat  sheet  stock  and  trim  to  the 
required  length  and  width. 

2.  Roll  form  each  strip  to  a  hat  shape. 

3.  Locate  the  uncured  hat  sections  over  a  rubber  mandrel. 

4.  Locate  the  rubber  mandrels  on  the  pre-cured  solid  laminate 
pieces  and  cure  in  place. 

5.  Remove  the  rubber  mandrels. 

The  concept  of  deriving  many  parts  from  a  single  master  does 
not  end  with  this  rather  simple  application.  GO  looks  forward  to 
manufacturing  wing  skins  in  almost  the  same  manner.  The  tape-laying 
machine  can  be  programmed  to  have  the  head  discharge  tape  only  at 
specified  intervals,  and  at  the  proper  orientation.  In  this  way. 


13 


complex  reinforcing  areas  can  be  built  up  on  a  flat  skin  and  then 
covered  with  the  outer  layers.  A  cutting  device  similar  to  a  cookie 
cutter  would  then  stamp  out  a  number  of  fabricated  skins  from  the 
master  lay-up.  As  stated  previously,  this  concept  has  not  been  used 
in  production  as  yet  since  it  is  predicated  on  the  availability  of 
effective  composite  cutting  tools. 

The  automated  tape  laying  machine  itself  is  far  from  perfected, 
however,  and  its  real  value  as  a  production  instrument  will  only  be 
realized  after  some  critical  problems  are  solved.  The  light-sensitive 
defect  detecting  apparatus  is  highly  sophisticated  and  doubts  have 
been  raised  as  to  its  effectiveness  in  spotting  a  defect  and  rejecting 
the  tape  when  the  machine  is  operating  rapidly.  The  variations 
from  supplier  to  supplier  in  the  tape  itself  make  it  difficult  to 
design  a  head  capable  of  accepting  every  brand  of  tape  available. 

The  amount  of  computer  programming  necessary  to  insure  that  the 
head  lays  the  correct  length  of  tape  at  the  proper  orientation 
within  specified  tolerances  is  a  factor  that  will  have  to  be 
reckoned  with.  Finally,  despite  the  fact  that  advanced  tape  laying 
heads  can  now  move  on  five  axes,  there  is  still  some  question  as  to 
their  ability  to  lay  up  pieces  with  a  high  degree  of  curvature  or 
complex  contours. 
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l.S.l.l  Hand  Lay-Up 

Because  of  the  above  problems,  GD  does  the  vast  majority  of  # 

its  tape  laying  by  hand.  The  material,  in  the  form  of  3"  tape, 
is  inspected  for  broken  and  crossed  fibers  and  gaps  between  the 
fibers  before  it  is  used.  The  tape  is  laid  up  on  a  metal  tool 
with  each  layer  at  the  proper  orientation.  Each  strip  of  tape 
in  the  same  layer  must  be  within  0.03"  of  the  adjoining  strips 
and  overlapping  is  not  permitted.  This  is  a  more  tedious  process 
than  using  a  machine  but  quality  control  as  better  and  complex 
curves  can  be  laid  up.  Most  of  GD's  work  has  been  experimental, 
such  as  the  F-5  fuselage  and  some  fuselage  bulkheads,  but  F-111 
doublers  have  been  made  in  some  quantity  using  this  process. 

2.3.2  Grumman  Azaoipacz  MsMiod 

The  only  actual  "on-line"  production  of  composite  aircraft 
components  is  being  done  by  Grumman  Aerospace  of  Bethpage,  N.  Y. 

(GA).  Under  contract  from  the  Navy,  they  are  producing  skins  for 

the  F14-A  horizontal  stabilizer  from  boron  epoxy  tape.  GA  terms 

the  operation  semi -automated  because,  although  the  work  is  done 

by  hand,  the  tape  lay-up  is  performed  on  a  device  which  they  developed 

and  call  the  "flintstone  machine"  (Figure  3).  Unlike  the  technique 

used  at  GD  where  a  bottom  layer  is  placed  on  a  tool  or  form  and  subsequent 
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MOVING  BELTS 


layers  are  built-up  one  atop  the  other,  Grumman  uses  the  flintstone 
machine  to  lay-up  each  individual  layer  separately.  The  outline 
of  each  layer  and  its  orientation  is  determined  by  a  computer  analysis 
of  the  loads  in  the  stabilizer.  This  outline  is  permanently  traced 
on  a  mylar  template  and  oriented  to  match  the  proper  orientation  of 
the  tape  for  that  layer.  For  example,  if  a  layer  of  45°  was  specified 
the  outline  of  that  layer  would  be  rotated  45°  on  the  mylar  template. 
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The  template  is  fastened  to  the  flintstone  machine  and  the  operator 
moves  the  belts  so  that  the  mylar  is  in  a  position  to  receive  the 
first  strip  of  tape.  The  composite  tape  is  unwound  from  a  reel 
by  hand  and  placed  on  the  template  over  a  portion  of  the  table  that 
is  backlighted.  Here,  the  operator  can  easily  inspect  the  tape 
for  broken  and  crossed  fibers  and  gaps  between  the  fibers.  If  the 
section  of  tape  is  defective,  it  is  discarded  and  another  section 
unwound  into  place.  If  the  tape  is  good,  it  is  pressed  onto  the 
mylar  template  and  trimmed  at  the  edges  to  exactly  fit  the  outline. 

The  operator  then  presses  a  button  and  the  machine  moves  the  template 
exactly  three  inches  so  that  the  next  section  of  tape  can  be  laid 
up  within  0.03"  of  the  previous  one.  This  operation  continues  until 
the  entire  outline  on  the  template  is  laid  up  with  tape.  The  mylar 
template  is  removed,  the  layer  of  tape  covered  with  plastic  to 
protect  it,  and  the  next  template  placed  on  the  flintstone 
machine.  When  all  the  layers  of  a  respective  skin  have  been  laid 
up,  the  skin  itself  is  then  built  up  on  a  metal  tool.  The  mylar 
template  has  a  series  of  holes  which  fit  over  pegs  on  the  tool  to 
insure  that  each  layer  has  the  correct  orientation.  When  complete, 
the  tape  is  pressed  into  the  layer  below  it  and  the  mylar  peeled 
back.  The  entire  skin  is  bagged  and  sent  to  the  autoclave  to  cure. 

The  end  product  is  indistinguishable  from  one  that  is  laid  up 
a  strip  at  a  time,  but  Grumman  claims  two  significant  advantages  using 
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their  method.  First,  the  fiintstone  machine  lends  itself  more  readily 
to  rapid  production  since  it  is  semi -automated  and  a  number  of  them 
can  be  used  at  once  to  lay  up  different  layers  of  the  same  skin. 
Second,  the  backlighting  on  the  area  where  the  tape  is  to  be  laid 
gives  visual  inspection  of  every  piece  of  tape  in  the  skin  thus 
providing  a  high  degree  of  quality  control.  Presently,  Grumman  is 
manufacturing  eight  skins  per  month  and  they  hope  to  double  that 
figure  by  the  end  of  1972. 

The  boron  skins  cover  an  aluminum  honeycomb  which  decreases 
in  density  moving  out  toward  the  edges.  Boron  is  mounted  on  the 
structure  by  titanium  fingers  shown  in  Figure  4.  In  this  way, 
all  drilling  can  be  done  through  the  titanium  rather  than  encountering 
the  difficulties  of  cutting  boron.  The  entire  stabilizer  structure 
is  outlined  in  Figure  5. 


Fig.  4:  Titanium  multi-step  splice  cross-sections 


Fig.  4A:  Edge  splice  (dimensions  in  inches) . 
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Fig.  4B:  Main  splice.  Adhesive  is  Metibond  329. 


Fiaure  5.  F-14  BORON  STABILIZER 
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2.3.4  ?fwbl2M  ofj  Ma44  ?n.oductlon 


Despite  Grumman 's  apparent  success  in  fabricating  the  F-14 
stabilizer  skins,  there  are  some  problems  which  must  be  overcome 
before  composites  can  be  used  on  a  mass  production  scale.  Probably 
the  greatest  single  impediment  to  successful  fabrication  with 
composite  materials  is  the  present  lack  of  design/manufacturing 
coordination.  The  people  working  on  advanced  ideas  which  incorporate 
graphite  or  boron  epoxy  must  communicate  with  those  involved  in 
manufacturing  to  a  greater  degree.  The  most  efficient  and  promising 
design  is  worth  nothing  if  the  structure  cannot  be  made  in  the  shop 
or  carries  a  prohibitive  price  tag.  An  understanding  of  the  special 
problems  involved  in  the  fabrication  of  composite  parts  will  help 
designers  direct  their  efforts  toward  a  component  that  is  easily 
producible.  The  design/fabrication  interface  must  be  fully  developed 
before  the  mass  production  of  composite  structures  is  economically 
feasible. 

It  is  possible  that  many  of  the  composite  parts  manufactured 
presently  are  "overdesigned".  The  use  of  a  computer  to  determine 
the  optimum  size  and  orientation  of  each  layer  in  a  piece  is  often 
expensive  and  results  in  a  complex  pattern  that  is  difficult  to 
assemble  and  insure  by  quality  control.  The  use  of  less  complicated 
orientations  of  layers  may  result  in  a  part  that  suffers  only  a 


21 


small  performance  penalty  and  is  much  more  cost  effective. 

One  of  the  largest  hardware  problems  faced  by  manufacturers 
is  the  inconsistency  in  the  amount  of  tack  on  the  tape  surfaces, 
even  in  batches  from  the  same  supplier.  Insufficient  tack  on  the 
surface  may  cause  the  layers  to  separate  and  produce  gaps  in 
the  component.  Quality  assurance  is  lowered  and  sophisticated 
testing  devices  are  necessary  to  determine  the  location  and  extent 
of  these  gaps.  On  the  other  hand,  excessive  tack  makes  the  tape 
difficult  to  handle  while  laying  up  by  hand  and  the  quality  of  the 
component  may  be  impaired  by  the  presence  of  overlaps  and  broken 
sections  of  tape. 

2.3.5  kttdmcvtz  Manuj^actuJUng  AppAoach&6 

The  entire  concept  of  lay-up  is  predicated  in  the  industry 
selection  of  3"  tape  and  wide  goods  as  the  form  of  the  raw  material. 
This  arrangement  may  be  necessary  to  fully  utilize  the  performance 
benefits  that  accrue  from  a  structure  of  many  unidirectional  layers, 
each  very  strong  in  its  respective  direction.  But  for  components  in 
which  the  strength/weight  ratio  is  not  so  critical,  the  use  of  a 
resin  matrix  impregnated  with  randomly  oriented,  short  graphite 
or  boron  fibers  should  be  investigated.  It  is  quite  possible  that 
some  of  the  current  manufacturing  techniques  employed  by  the  fiberglas 
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industry  could  be  successfully  applied  to  composite  technology. 

The  most  promising  appears  to  be  the  sheet  molding  compound  process 
(SMC)  which  is  used  to  manufacture  automobile  parts  and  panels  such 
as  the  Corvette  body.  The  raw  material  is  made  by  chopping  long 
strands  of  glass  fiber  as  they  are  pulled  over  the  resin  matrix  so 
that  they  drop  onto  the  resin  in  a  random  manner.  The  compound  is 
rolled  onto  drums  and  ready  to  be  used.  Although  there  are  no 
time-temperature  requirements  such  as  those  mentioned  on  page  5 
with  fiberglas  in  this  form,  it  is  not  clear  whether  boron  or 
graphite  would  also  be  relieved  of  this  problem.  Parts  are  fabricated 
from  SMC  simply  by  unrolling  the  desired  length,  placing  it  on  a 
mold,  and  curing  it  by  compression.  Rapid  production  of  parts 
on  a  large  scale  using  this  method  is  economically  practical  as 
evidenced  by  the  number  of  fiberglas  firms  using  it.  Even  if  a 
structure  made  from  boron  or  graphite  in  this  form  could  not  meet 
performance  criteria,  SMC  could  still  prove  to  be  a  feasible  method 
for  utilizing  scrap  to  fabricate  pieces  which  are  not  critical  in 
carrying  loads.  From  a  fabrication  point-of-view,  an  investigation 
is  necessary  to  study  the  economic  implications  of  continuing  to 
manufacture  boron  or  graphite  parts  from  3"  tape. 


2.4  Curing  Process 


Once  the  part  is  laid  up  it  is  sent  to  the  autoclave  to  cure. 
Curing  permits  the  resin  to  redistribute  around  the  boron  or  graphite 
fibers  so  that  the  part  is  an  integrated  whole  rather  than  a  number 
of  layers  stuck  together,  A  typical  autoclave  setup  is  shown  in 
Figure  6. 


AUTOCLAVE  SHELL 


The  part  is  not  placed  directly  on  a  metal  tool  but  rather  on 
a  layered  bleeder  system  which  separates  the  two.  The  bleeder  system 
completely  encloses  the  composite  part  to  be  cured  and  is  necessary 
to  remove  the  excess  resin  which  flows  from  the  piece  at  elevated 
temperatures.  A  typical  bleeder  system  consists  of  the  following 
materials  starting  at  the  part  surface: 

1.  A  layer  of  TX1040  separator  to  prevent  the  bleeder  cloth 
from  sticking  to  the  part. 

2.  Layers  (depending  on  the  part  thickness)  of  paper 
bleeder  or  120  dry  glass  to  absorb  the  excess  resin. 

3.  A  layer  of  mylar  film  with  perforations  at  specified 
intervals. 

4.  A  layer  of  ply  paper  vent  cloth  to  measure  the  amount 
of  resin  coming  through  the  mylar  film.  This  is  a 
check  to  determine  if  the  part  is  being  bled  correctly. 

5.  A  plastic  pressure  bag  to  cover  the  entire  composite  part. 

One  of  a  number  of  different  cure  cycles  is  followed  depending 

on  the  type  of  material  used.  For  example,  the  basic  cure  cycle 
for  Narmco  5505  specifies  2  hours  in  the  autoclave  at  85  psi  and 
350°F.  However,  minor  changes  are  frequently  made  to  compensate 
for  part  thickness,  resin  flow,  and  few  other  small  details. 
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There  are  a  number  of  problems  involved  in  the  curing  process 
which  are  presently  being  examined.  The  most  significant  hardware 
problem  is  the  volume  content  and  type  of  resin  used  in  the  matrix. 

An  incorrect  estimate  of  the  resin  volume  or  chemical  composition 
can  cause  serious  problems  in  bleeding  during  cure.  The  resin  may 
flow  too  quickly,  too  slowly,  unevenly,  or  possibly  not  at  all, 
all  of  which  can  cause  the  part  to  be  defective. 

One  problem  in  the  cure  cycle  that  would  probably  not  be 
readily  forseeable  is  the  tedious  task  involved  in  placing  and 
removing  the  plastic  vacuum  bag  over  the  part.  Placing  the  tape 
around  the  edges  is  time-consuming  in  itself,  but  removal  of 
the  tape  and  bag  after  curing  has  taken  place  is  even  more  difficult. 
Toward  eliminating  this  problem,  Grumman  is  doing  development  work 
on  improved  advanced  composites  bag  molding  processes.  They  have 
developed  a  silicone  blanket  with  a  vacuum  actuated  0-ring  to 
hold  the  blanket  to  the  part.  Work  is  not  complete  on  the  project 
but  the  major  advantages  that  appear  feasible  using  the  silicone 
blankets  are  reduced  lay-up  and  removal  time,  re-usable  bagging, 
and  applicability  to  complex  shapes.  Successful  development  of 
such  a  process  would  be  a  significant  step  toward  the  fabrication 
of  composite  parts  on  a  mass-production  scale. 
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The  greatest  expense  involved  in  the  curing  process  is  the 
time  required  for  the  cure  cycles  while  the  part  is  in  the  autoclave. 
Grumman  is  currently  conducting  a  test  program  using  microwaves  at 
the  frequency  of  FM  broadcasting  (-100  MHz);  this  could  cut  cure 
time  and  autoclave  expense  tremendously.  Tests  have  been  performed 
on  fiberglass  reinforced  resins  with  satisfactory  cure  being  achieved 
in  less  than  half  an  hour.  When  this  cure  time  is  compared  to  those 
involved  with  heat-treating  processes  the  time  saved  can  be  measured 
in  houM.  Grumman  is  presently  seeking  additional  funding  from  the 
Air  Force  Materials  Lab  to  run  a  full  scale  testing  program  on 
boron  and  graphite  impregnated  resins.  It  is  quite  obvious  that 
if  boron  or  graphite  epoxy  can  be  cured  using  microwave  techniques, 
the  cost  of  the  finished  component  could  be  significantly  reduced. 
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2.5  Machining  Problems 


A  serious  drawback  associated  with  the  use  of  advanced  composites 
is  the  difficulty  involved  in  machining  the  finished  product.  Parts 
are  usually  laid  up  to  perfectly  specified  dimensions  so  that 
trimming  of  the  composite  material  itself  is  not  necessary. 

However,  there  are  locations  where  bolt  holes  and  other 
mounting  points  must  be  drilled.  Because  of  the  hardness  of  the 
fibers  in  advanced  boron  composites,  conventional  drilling  and 
machining  methods  deliver  only  minimal  results.  Diamond- tipped 
drills  have  shown  some  degree  of  success  but  these  tend  to  wear 
out  relatively  quickly. 

Some  manufacturers  have  taken  steps  in  the  fabrication  process 
to  eliminate  the  necessity  of  drilling  a  composite  skin.  Grumman 
uses  the  stepped  titanium  fingers  shown  in  Figure  4  to  attach  the 
boron  skin  to  the  stabilizer  substructure  because  mounting  points 
can  be  drilled  more  easily  in  the  titanium.  There  has  been  discussion 
about  the  prospect  of  laying  up  the  composite  around  plugs  which 
will  be  removed  after  cure  to  provide  the  proper  holes.  The 
possibility  also  exists  that  the  plugs  will  be  metal  which  can  be 
drilled  out  after  cure  to  provide  mounting  points  similar  to  the 
titanium  finger  concept  utilized  by  Grumman. 
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The  real  breakthrough  in  advanced  composite  machining,  however, 
appears  to  be  in  the  area  of  ultrasonic  drilling.  Funded  by  the 
Air  Force,  Grumman  is  currently  experimenting  with  ultrasonic 
techniques  on  boron/epoxy.  Thus  far,  the  results  have  shown 
excellent  potential  for  the  use  of  ultrasonic  equipment  to  signifi¬ 
cantly  reduce  the  timer? and  difficulty  involved  in  drilling  advanced 
composites.  A  special  diamond-studded  drilling  bit  has  demonstrated 
good  endurance  when  used  with  ultrasonic  assistance.  Since  most 
of  the  structures  being  fabricated  are  composite  skin  and  honeycomb 
metal  core  it  is  of  particular  importance  that  drill  bits  which 
work  well  on  metal  are  poor  for  advanced  composites  and  vice  versa. 
When  drilling  a  simple  hole  through  such  a  structure,  the  metal- 
favoring  bit  is  used  and  the  penalty  incurred  on  the  thin  skins  is 
accepted.  However,  if  a  countersunk  hole  is  necessary,  using  only 
one  type  of  drill  bit  is  unacceptable.  To  alleviate  this  problem, 
Grumman  has  developed  the  two-stage  drill  bit  shown  in  Figure  7. 

The  neck  is  a  metal -favoring  composition  which  is  hindered  as  it 
drills  initially  through  the  composite  skin  but  becomes  very 
effective  in  the  metal  core.  The  tapered  base  is  a  composite- 
favoring  composition  which  allows  the  countersink  to  be  drilled 
in  the  skin  in  one  operation.  This  combination  drill  bit  coupled 
with  ultrasonic  assistance  promises  to  significantly  reduce  some 
of  the  machining  time  involved  in  the  fabrication  process. 


29 


II  I 


CMC  M  i  C  A  1. 


6r<)«ichev  e«.-*oft.  M.-jis  •  Donburv,  Contt. 

Pkilcitii'littMo  Pp  •  SfIvcJ  Spring!, 


Figure  7.  COMBINATION  DRILLING  -  COUNTERSINKING  TOOLS 
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2.6  Conclusions 

The  Air  Force  and  all  the  various  advanced  composite  manufacturers 
have  invested  a  tremendous  amount  of  time  and  money  in  amassing  the 
current  store  of  knowledge  on  composite  materials.  The  fact  that 
they  are  now  becoming  concerned  with  the  direct  problems  of  fabrication 
demonstrates  that  the  materials  are  no  longer  considered  just  an 
experimental  venture.  Prior  work  and  investigation  have  developed 
a  confidence  In  composite  materials  as  an  effective  substitute 
for  conventional  metals  in  aircraft  construction.  The  focus  has 
now  been  placed  on  the  practicality  of  advanced  composites  on  a 
cost  basis  rather  than  a  performance  or  reliability  one.  There  are 
many  problems,  some  of  them  mentioned  in  the  report,  which  remain 
unsolved  in  the  realm  of  composite  fabrication.  But  the  fact  that 
companies  are  addressing  their  efforts  toward  new  tooling  materials, 
automated  fabrication  techniques,  and  sophisticated  yet  cost  effective 
practices  in  curing  and  machining,  clearly  demonstrates  that  advanced 
composite  materials  have  a  tremendous  potential  in  the  iirmediate 
future. 
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CHAPTER  III 


A  COMPUTER  PROGRAM  FOR  WEIBULL  STATISTICS 
3.1  Review  of  Weibull  Statistics 

The  following  report  outlines  the  basic  equations  of 
Weibull  statistics  and  presents  a  computer  program  developed 
in  order  to  fit  empirical  data  to  the  extreme  value  (Weibull) 
distribution.  The  computer  program  takes  unsorted  numerical 
data  and  performs  all  of  the  standard  sorting  and  cataloging 
routines  necessary  for  statistical  analysis.  The  data  is  then 
fit  in  a  least  square  sense  to  a  Weibull  distribution  function 
and  a  chi-square  confidence  test  is  employed. 

The  basis  for  a  fatigue  life  prediction  includes 
statistical  inference  from  an  empirical  data  base,  and  theoretical 
or  engineering  reasoning  to  choose  a  probability  model.  It  is 
now  recognized  that  the  extreme  value  (Weibull)  distribution  is 
the  best  physical /theoretical  model  for  the  fatigue  process. 

Generally  speaking,  the  shape  parameter  (scatter)  may  be  established 
by  component  testing  and  the  scale  parameter  (roughly,  the  mean) 
may  be  established  by  a  few  (1-3)  full  scale  tests. 

The  analytical  reliability  (or  1 -cumulative  distribution) 
cumulative  function  for  the  Weibull  distribution  is  given  by 

R(x)  =  exp  -  [(x/p)"^]  (1) 

where  a  is  the  shape  parameter  and  p  is  the  scale  parameter.  The 
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mean  value  for  any  distribution  density  function  f(x)  is  given  by 

00 


X  = 


xf(x)dx 


and  the  variance,  by 


r2  = 


=  J*  (x-x)2  f(x)dx 

—00 

CO 

=  /x^f(x)dx-x 

—00 

For  the  Wei  bull  distribution  (1),  the  mean  is  given  by 

X  =  e[r(l  +  1/a)] 


and  the  variance  is  given  by 

a2  =  e2[r(l  +  2/a)  -  r2(l  +  l/a)] 

In  general,  the  coefficient  of  variation  is  a/x  such  that  for  (1) 

cv  =  a/x  = /TTn~+~27^  r^(l  +  l/a)J  /  r(l  +  1/a) 

In  treating  experimental  results  one  chooses  plotting 
position  reliability  (e.g.  see  Freudenthal  [10]  Gumbel  [11]). 

For  the  m^^  ordered  data  point  in  a  sample  of  size-n,  the  empirical 
reliability  is  given 

R(m)  =  1  -  m  /  (n  +  1 ) 

The  median  for  (7)  is  found  by  plotting  R(m)  vs.  x^^  and  taking  the 


intercept  at  R  =  0.5.  The  variance  from  the  empirical  data  is 
computed  from  the  mean  (average) 


The  median  time  to  first  failure  is  denoted  x,  in  (1), 
and  in  a  finite  sample  by  m  =  1  in  (7)  such  that 

exp  -  [(xi/3)“]  =  1  -  1  /  (n  +  1) 

which  can  be  solved  for  x^ 

Xi  =  3  [-£n(l  -  ^ 

It  has  been  found  from  extreme  value  statistics  that  the 
distribution  of  first  failures  in  samples  of  equal  size-n  is  the 
same  as  the  parent  population;  i.e.  the  same  functional  form 
and  the  same  a.  This  is  illustrated  for  the  normal  distribution 
in  Figure  8. 


1  i  f  e  ti  me 


Figure  8:  First  Failure  Distributions 


The  scale  and  shape  parameters  for  the  first  failure  distribution 
are  given  by 

3l  =  Bn‘‘'/“  ;  ai  =  a  (11) 

The  frequency  density  function  for  first  failures  in  samples 
of  size-n  is  given  by 

fi  (x)  =  n  f(x)  exp  (V')(i7)  ] 

Eq.  (12)  can  be  derived  on  the  following  basis.  Let  R(x)  be  the 
probability  of  survival,  i.e.  l-F(x),  of  a  single  item.  The 
probability  of  survival  of  n-i terns  is  [R(x)]'^,  v/hich  is  the 
probability  of  no  failures  in  n-i terns.  The  probability  of  failure 
of  one  out  of  n-i terns  is  then 

Fi(x)  =  1  -  [1  -  F{x)]  "  (13) 

from  which  the  probability  density  function  is  obtained 

fi(x)  =  n(l  -  F(x))""''  f(x)  (14) 

Eq.  (14)  for  the  Wei  bull  density  gives  (12). 

The  following  sections  outline  the  operation  of  the 
computer  program  as  well  as  the  numerical  and  mathematical 
foundations  for  each  section.  The  computer  program  listing  in 
a  sample  problem  output  are  included. 
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3.2  Detailed  Program  Description 
3.2.7  kfvtcuj  VlmzYi&loyHi 

All  arrays  are  given  a  dimension  of  100.  This  is  arbitrary, 
and  can  be  changed  by  altering  cards  (1)*,  (2),  and  (3).  In  general,  the 
following  rules  must  hold  pertaining  to  array  dimensions.  The  following 
arrays  must  have  dimension  greater  than  or  equal  to  the  number  of  data 
points;  X,  Z,  F,  FI,  Fll,  and  P.  The  arrays  B1 ,  B2,  VM,  T,  and  T1 
must  have  dimension  greater  than  or  equal  to  the  number  of  groups 
desi gnated. 

3.2.2  Vdta  Po-inti>,  Gfioiip  NimbeJti>  (1)  [12] 

The  number  of  data  points  is  input.  This  includes  the  total 
number  of  data  values,  even  though  two  or  more  values  may  be  the  same. 

The  data  values  may  be  input  in  any  convenient  order.  The  number  of 
groups  is  selected  for  the  purpose  of  setting  up  a  histogram.  That 
is,  any  convenient  number  of  groups  may  be  chosen.  In  general,  to 
make  a  histogram  meaningful  a  large  number  of  data  points  is  necessary. 
The  number  of  groups  would  increase  with  the  number  of  data  points  with 
a  maximum  of  perhaps  ten.  The  correlation  between  number  of  groups 
and  number  of  data  points  is  quite  arbitrary. 

*Parentheses  refer  to  program  operation  sequence  numbers  as  listed  in 
Section  3.4. 
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The  input  formats  for  N  (number  of  data  points)  and  N2  (number 
of  groups)  are  identical.  They  consist  of  printing  three  digits  on  the 
first  3  positions  on  two  cards.  Both  must  be  right  justified  and  are 
integers.  No  decimal  point  is  necessary. 

Data  values  are  input  as  real  numbers.  They  are  also  right 
justified  with  decimal  points.  Each  value  is  allotted  ten  positions 
for  digits  and  decimal  point.  A  maximum  of  five  values  is  allowed  on 
each  card. 

The  selected  critical  chi-square  value  (generally  used  for 
comparison  with  the  calculated  value  and  designated  by  CHCRIT)  and 
its  associated  confidence  level  (SI6LEV)  are  input  identically.  They 
are  on  two  separate  cards,  real  numbers,  and  are  right  justified  with 
decimal  points.  The  first  ten  spaces  on  these  two  cards  are  allotted 
for  these  values.  The  CHCRIT  value  and  its  associated  SIGLEV  are 
chosen  with  reference  to  the  number  of  data  points.  For  N  equal  to 
the  number  of  data  points,  N-1  is  equal  to  the  number  of  degrees 
of  freedom  of  the  test  group.  Generally  tables  of  values  of  N-1  vs. 
SIGLEV  result  in  values  of  CHCRIT.  These  tables  are  found  in  many 
mathematical  reference  books. 

In  summary,  input  values  follow  the  following  form: 


CARP 

INPUT 

FORM 

Card 

1 

No.  of  data  values  (N) 

(13) 

Card 

2 

No.  of  groups  (N2) 

(13) 

Cards 

►  2r^ 

data  values  (X(I)) 

5F10.0 

Card 

4+ 

selected  crit  (CHCRIT) 

FIO.O 

Card 

5+ 

selected  conf  lev  (SIGLEV) 

FIO.O 
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In  the  case  of  N  and  N2,  the  maximum  value  inputted  is  999.  To  increase 
this,  the  format  must  be  changed. 

For  all  input  parameters  with  formats  of  FI 0.0,  the  maximum 
input  value  is  999,999,999.  To  increase  this,  the  format  must  be 
changed.  Data  values  are  represented  in  the  array  X. 

3.2.3  Seizing  Data  VaJiau 

The  data  values  are  ranked  in  ascending  order,  i.e.  X(l)  = 
smallest  value,  X(N)  =  largest  value.  This  is  done  to  simpligy 
operations  carried  out  later  which  require  differentiation  of  the 
relative  magnitude  of  the  values.  Cumulative  frequency  evaluation, 
for  example,  requires  that  its  value  (F(l))  grows  as  the  data  value 
grows . 

The  method  of  sorting  is  known  as  a  'bubble  sort'  where  the 
following  procedure  is  used.  The  first  two  data  values  are  compared 
and  the  largest  chosen.  It  is  then  compared  to  the  next  value  and 
the  largest  of  those  two  is  chosen.  It  is  then  compared  to  the  third,  etc. 
Finally  the  largest  value  remains.  The  procedure  is  repeated  for  I  =  1 
through  N-1  values  this  time.  The  remaining  value  will  now  be  the 
second  largest.  The  procedure  is  repeated  until  all  values  are  ranked. 
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3.2.4  CaZaatation  Cumutoutivz  Tfitqiimcij  (4)  [13] 

The  cumulative  frequency  value  for  each  data  point  is  calculated. 
For  the  present  study  the  Wei  bull  distribution  will  be  used  to  describe  both 
the  frequency  density  and  cumulative  frequency  functions,  as  given  by 


cumulative  frequency 

F(x)  =  1  -  e  "(?) 

frequency  density 


(15) 


f(x) 


dF(x)  _  _a  ot-l 


(0 


dx  6“ 

where  a  and  3  are  the  shape  and  scale  parameters  respectively.  It  should  be 


noted  that  these  functions  are  valid  for  x  ^  0,  a  characteristic  which 
determines  lower  bounds  on  groups  chosen  for  a  histogram  study. 

The  cumulative  frequency  for  each  data  point  is  given  by  the 
following  expression 

F(X. )  =  i/T+N  (16) 

where  F  is  the  cumulative  frequency  associated  with  X^. ,  the  ith  data 
point  in  ascending  order.  N  is  the  total  number  of  data  points. 

This  particular  choice  of  F(X)  forces  the  function  to  pass 
through  the  median  which  forces  an  equal  number  of  data  points  to  fall 
above  and  below  the  median  value  on  a  graph  of  F(X)  vs.  X.  Cumulative 
frequency  values  are  represented  in  the  array  F. 
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3.2.5  CaZcutatlon  OjJ  Statutlc-oZ  Pcuiam^eJiA  (6-8,  10,  11)  [12,14] 

Several  values  which  are  of  interest  when  studying  a  particular 
distribution  are  the  range,  mean,  median,  standard  deviation,  and  co¬ 
efficient  of  variation.  The  range  is  calculated  by  taking  the  diff¬ 
erence  of  the  largest  and  smallest  data  vlaues.  That  is,  R  =  X(N)  - 
X(l)  for  ascending  values.  The  mean  or  average  is  the  sum  of  the 
data  values  divided  by  the  number  of  data  points.  It  is  represented 
by  Q. 

The  median  of  a  ranked  group  of  numbers  is  defined  as  the 
middle  value  or  the  mean  of  the  two  middle  values.  It  is  in  a  position 
such  that  the  number  of  data  values,  which  are  greater  than  its  value 
is  equal  to  the  number  of  data  values  which  are  less  than  its  value. 

For  N  odd,  the  median  would  be  given  by  X(N+l/2).  For  N  even,  the 
median  is  given  by  the  average  of  X(N/2)  and  X(N/2+l).  It  is  rep¬ 
resented  by  DMED  in  the  program. 

The  standard  deviation  (STDDEV)  is  given  by: 

N 

a  =  (X.-X 

i-1 

where  X.  is  the  ith  data  value  and  X  is  the  mean. 

1  m 
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3. 2.  (5  Gao  up  (fJldtk  Calcutcutcon 

The  group  width  designates  the  width  of  each  section  used  on  a 
histogram  plot.  Its  size  is  arbitrary  and  the  following  discussion 
outlines  the  reasons  for  its  particular  selection  in  this  program. 

AX,  the  width,  designated  by  STEP  in  the  program  is  simply 
given  by  X(N)/N2,  the  largest  data  value  divided  by  the  desired  number 
of  groups.  Another  method  of  selection  might  be  defining  aX  as  the 
range  divided  by  the  desired  number  of  groups.  A  histogram  is  a  very 
arbitrary  plot  in  any  case,  and  the  width  chosen  here  seems  to  be  as 
appropriate  as  any  other. 

3.2.7  Gaoup  BoundcuvieJi  CalauLatO-d;  NambeA  ol  VaZa.  FolyvU  -in  Each.  GAoup 
VeXoAmns-d  {14,  15) 
liid-HalueA  ol  Gao  ups  CalcuZattd  ilS] 

Selection  of  group  boundaries  for  a  histogram  plot  is,  like 
the  group  width,  a  somewhat  arbitrary  choice.  That  is,  as  far  as 
selecting  the  lower  boundary  limit  of  the  first  group  is  concerned. 

After  that  selection  is  made  all  other  boundaries  are  defined  by  the 
group  width. 

The  lower  limit  of  the  first  group  is  taken  to  be  0.0.  This 
means  that  all  data  values  must  be  greater  than  0.0  which  is  a  character¬ 
istic  of  the  Weibull  distribution.  Other  statistical  functions  might 
have  data  values  lying  on  both  sides  of  the  dependent  variable  axis. 

The  value  of  STEP  is  added  to  0.0  to  get  the  first  upper  boundary,  and 
each  succeeding  limit  is  determined  by  adding  the  value  of  STEP  again. 
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Each  data  point  is  then  checked  to  see  which  group  it  lies  in 
and  the  total  number  of  data  points  in  each  group  is  recorded.  In  some 
instances  computer  round-off  errors  cause  less  data  points  to  be  included 
in  the  last  group  than  actually  exists.  To  check  for  this  the  total 
data  points  in  the  groups  is  determined  by  summing  the  number  of  points 
in  each  group.  If  this  sum  is  less  than  N,  the  total  number  of  data 
points,  enough  points  are  added  to  the  last  group  sub-total  to  correct 
the  error. 

When  checking  to  see  whether  or  not  a  data  value  falls  between 
two  boundaries  the  upper  limit  is  inclusive  whereas  the  lower  limit  is 
not  (except  for  0.0,  which  is  inclusive).  The  mid- value  of  each  group 
is  found  by  averaging  the  boundaries  of  the  group.  The  number  of  data 
points  in  each  group  is  represented  by  the  array  T,  whereas  mid-values 
are  represented  by  VM. 

3.2.  S  Shape,  and  Scate.  PaJiameteAA  Calculatzd  (16,17) 

In  order  to  solve  for  a  and  3,  the  shape  and  scale  parameters, 
the  following  relationship  is  used: 

In  {  In  {1/(1  -  F(X))  }}  =  alnx  -  aln  3  (18) 

which  is  easily  derivable  from  the  defining  equation  for  F(x).  In  this 
expression  X  represents  an  individual  data  value  and  F(X)  is  its  associated 
empirical  cumulative  frequency  value.  That  is,  for  X=X^. ,  F(X)  =  i/N+1. 

This  relationship  indicates  that  In  {  In  {1/(1-F(X))  }}  is  a 
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linear  function  of  InX.  Therefore  if  the  two  functions  are  plotted  on 
an  appropriately  scaled  graph  a  straight  line  will  result.  In  order 
to  determine  a  and  3  a  least-squares  curve  fit  will  be  used. 

The  first  step  is  to  transform  F(X)  and  X  into  their  respective 
functions  of  the  natural  logarithm.  Inspecting  the  relationship  again, 
it  can  be  shown  that  the  slope  of  the  line  is  given  by  a  and  the  inter¬ 
cept  is  -a£n3.  Let 

F^.  =  £n  {  £n  {  1/(1-F{X.))}}  (19) 

F2-  =£n(X.)  (20) 

Then  for  y  =  mx  +  b,  y  =  F2^ ,  X  =  F-j^. ,  m  =  a,  and  b  =  -aln^.  The  follow¬ 
ing  relationships  are  true  for  the  linear  curve  fit  procedure: 


■oi£n3 


N  N  N  N  p  N  2  , 

=  [.I:(F2P'  E  F2,F„]  /  -  (EF^P']  (2.9) 


If  the  righthandside  of  the  above  expression  is  given  as  Q,  then 

3=  l/EXP(Q/a)  (21) 

3.2.9  Tlieo/Lc££aa£  F/Leqaency  o{j  Individual.  Vaia  ?ointi>  (19) 

The  individual  frequencies  for  each  data  point  are  derived  from 
the  theoretical  distribution.  What  this  value  represents  for  each  point 
is  the  number  of  times  it  would  appear  in  a  group  of  N  specimens  as  pre¬ 
dicted  by  the  Weibull  distribution  defined  by  the  a  and  3  previously 
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calculated.  These  frequency  values  are  known  as  the  expected  frequency 
values  and  will  be  used  in  a  goodness-of-fi t  calculation  later  in  the 
program. 

Consider  the  frequency  density  function,  f(x).  It  has  been 
defined  in  one  way  as 

f(x)  =  dF(x)/dx  (22) 

The  upper  limit  of  this  function  is  ^  and  the  lower  limit  0.  The 
following  is  therefore  true: 

CO 

/f(x)  dx  =  F(c»)  -  F(0)  =  1.0  (23) 

0 

That  is,  the  total  cumulative  frequency  is  100%  or  1.0.  This 
total  is  made  up,  theoretically,  of  a  summation  of  an  infinite  number 
of  f(x)  values  multiplied  by  an  infinitesimal  width  dx.  Actually,  we  have 
a  finite  number  of  f(x)  values  (that  is,  a  finite  number  of  X  values). 
Because  of  this,  it  follows  that  dx  becomes  aX,  a  finite  width.  Therefore 
there  actually  exists  a  series  of  f(X^. )  AX^.  v/hich,  when  summed,  total 
1.0,  the  cumulative  frequency  over  the  entire  range  of  X^  values. 
Theoretically 

b 

f(x)dx  =  F(b)  "  F(a)  (24) 

a 

This  is  illustrated  by  the  graph  given  below: 
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f(x) 


of  data  values  which  falls  between  Xj^  and  x^.  That  is,  given  an  infinite 

population  the  percentage  of  values  which  fall  between  X|^  and  x^  is  equal 

to  F(b)-F(a).  For  a  finite  population,  this  area  will  be  used  to  represent 

the  percentage  for  a  single  data  value.  For  the  above  graph  this  value 

would  be  X..  The  question  is:  "For  an  X.  in  a  population  of  N  values, 

how  are  x  and  x.  determined  so  as  to  give  a  good  representation  of 
a  b 

F(b)-F(a)?"  The  choice  is  quite  arbitrary,  and  the  following  has  been 
chosen  for  use  in  this  program.  Consider  the  graph  below: 


Figure  10:  Frequency  Groups 
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For  any  x. ,  a  given  data  value,  and  Xj^,  the  limits  which  determine 
F(a)  and  F(b),  are  given  as  the  mid-points  between  x^  and  x^_i,  and 
x^.^1  and  x^  respectively.  This  would  imply  that 

„„.p(a,=  (35, 

=  expected  frequency  of  x^. 

_  )<  “ 

For  F(x)  =  1  -  e  e 

F(b)-F(a)  -  e  "  (^)  -  e  “(^) 

This  is  the  calculation  method  used  by  the  program  to  give 
expected  frequency  values  for  individual  data  points.  The 
individual  frequencies  are  respesented  by  the  array  P  and  are 
listed  under  'CALC  FREQ'. 

3.2.10  Ccitc.atatQ,d  V^^qumcy  Each  GAoup  (20) 

These  values  represent  the  total  percentiles  for  the  groups 

defined  previously.  All  they  represent  are  the  summation  of  the  individual 

frequencies  calculated  over  the  total  data  points  in  any 

particular  group. 

The  program  checks  each  group  to  see  if  it  contains  data  points. 
The  first  group  to  contain  data  points  may  have  m  values,  m  ^  n.  Then 
the  first  m  individual  frequencies  are  summed.  The  next  group  to  have 
a  non-zero  number  of  data  values,  say  £  values,  then  has  the  next  I 
individual  frequencies  summed. 
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These  values  are  represented  by  Fll(I),  and  are  printed 
out  in  'Calc  Freq  in  Group.' 

3.2.11  Ob^eJLVud  FA.zqamcy  Catculation  Each  Gwap  (21) 

For  each  data  point,  the  observed  frequency  is  given  by 

1/N+l.  This  is  consistent  with  the  choice  of  F(X^.)  =  i/N+1  for 
cumulative  frequency  value «  For  a  group,  the  total  number  of 
data  values  in  that  group  is  determined  and  multiplied  by 
1/N+l.  These  values  are  represented  by  T1(I)  and  are  listed 
under  'OBS  Freq  in  Group'. 

3.2.12  Chl-SqaaJifL  CatcutatLon  (23)  [14] 


The  goodness-of-fit  test  which  will  be  used  in  this  program 
is  the  Chi-Square  (x^)  test.  Although  for  a  smaller  number  of  data 
points  the  Komarov-Smi rnoff  test  may  be  better,  the  Chi-Square  test 
has  been  chosen  due  to  its  wide  applicability  and  the  relative  ease 
by  which  the  value  of  Chi-square  and  its  associated  confidence  level 


can  be  calculated. 


The  value  of  tor  a  given  function  is  defined  as: 

...  ^  (27) 


i=l  '-i 

where  0^.  =  observed  frequency  for  data  point  i  (=1/N+1) 

e^.  =  expected  frequency  for  data  point  i 

This  value  is  represented  by  CHISQ  and  is  printed  out  under  the  title 
'Calc  Chi-Square' . 
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3.2.13  Coni^^dence.  Leuet  CalcuZated  Chl-SquoA^  (24) 


In  order  to  determine  the  confidence  level  for  which  the 
calculated  chi-square  is  a  critical  value  the  area  beneath  a  chi-square 
curve  is  determined.  If  the  total  area  under  the  curve  is  normalized  then 
the  area  under  any  portion  from  0  to  some  value  of  represents  the 
confidence  level  of  any  value  of  and  an  associated  number  of  data 
points  is  given  by  the  following  expressions; 

For  N  even: 


Q(x^|N)  =  2Q{x)  +  sZ(x) 


yZr-l 

l*3-5...(2r-l) 


(28) 


X 


For  N  odd: 


Q(x^iN)  =  ^  Z(x) 


N-3 

2  2^ 

S  2.4^.V(2r) 
r=l  '  ' 


(29) 


where:  Q(x^l!^)  represents  the  confidence  level 


Z(x) 


_1_  e 


Q(x) 


j  Z(t)dt 
X 


(30) 


(31) 


Because  the  form  of  Q(x)  given  above  is  in  terms  of  a  non-integrable 
function,  the  following  approximation  is  made: 
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Q(x)  =  Z(x) 


I_  1_  1_  i- 

x+  x+  x+  x+  x+ 


9 


(32) 


a  continued  fraction  expansion.  This  is  expressed  as  the  following 


1 

X  +  1 

X  +  2 

X  +  3 


(33) 


This  infinite  summation  is  used  to  calculate  Q(  )  in  the  program.  At 
this  time  the  error  limit  is  given  as  .0001.  If  this  is  found  to  be 
too  large  it  can  easily  be  altered  (card  232). 

In  summary  then,  the  following  is  used  to  calculate  the 
confidence  levels: 

N  even: 


Q(xVN)  =  - 

I_  3_  4_  ^  \ 

X+  X+  X+  X+  / 

N-2 

2r-l 

(34) 

+  2  -xV2 

L 

E 

r=l 

X 

-I-  — —  e 

l*3-5...(2r-l) 

Q(xVN)-  (35) 

2 

Q(x  /N)  is  represented  by  Q2  in  the  program. 

3.2,14  Compcuilion  CaZaalat^d  and  SdLzcXtd  Chd-Squxuit  UoZiiu  (26) 

The  calculated. value  of  chi-square  is  compared  to  the  selected 
critical  value  of  chi-square.  This  can  also  be  done  by  the  operator  as 
both  values  are  printed  out  along  with  the  calculated  confidence  level. 
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After  comparing  the  two  numbers,  a  conment  is  printed  out  regarding  the 
comparison.  Basically  the  comments  read  as  follows: 

1)  if  the  calculated  chi-square  is  the  larger  of  the  two 
then  at  the  selected  confidence  level  for  the  N-1  degrees 
of  freedom  the  distribution  found  by  the  program  is  not 
good  enough.  This  probably  means  that 

a)  too  few  data  points  exist  for  a  good  fit  or 

b)  the  scatter  is  such  that  the  Weibull  distribution 
is  not  an  appropriate  descriptive  function 

2)  if  the  calculated  chi-square  is  the  smaller  of  the  two 
than  at  the  selected  confidence  level  for  the  N-1  degrees 
of  freedom  the  distribution  found  by  the  program  is  a  good 
description.  (One  might  note,  however,  that  a  calculated 
chi-square  which  is  too  small,  and  consequently  a  confidence 
level  about  1.0,  may  also  indicate  something  is  wrong  with 
either  the  distribution  or  the  data  selection.) 


50 


3.3  Program  Description 


In  order  to  give  the  reader  a  quick  but  clear  description  of 
the  attached  computer  program,  the  following  will  be  used  as  a  presentation 
guideline: 

1)  The  program  operations,  either  I/O  or  internal  functions, 
will  be  listed  in  the  order  in  which  they  are  performed. 

2)  These  operations  will  be  accompanied  by  card  numbers 
corresponding  to  the  program  listing. 

3)  A  complete  program  description  will  be  given.  The  major 
operations  on  the  sequence  list  will  be  described  as  to 
how  the  operation  is  done  and  why  it  is  done  that  way. 

In  the  case  of  input  the  format  used  will  be  described 
in  detail . 

4)  Immediately  following  the  description  titles  will  be 
given  the  item  number(s)  of  the  appropriate  sequence 
listing. 

5)  Following  item  number(s)  will  be  a  list  of  possible 
references  for  further  discussion.  (Each  description 
will  not  necessarily  have  a  reference  listed.)  The 
references  are  given  by  numbers  and  listed  together 
at  the  end. 

In  general,  each  major  operation  or  group  of  similar  operations 
(internal)  will  be  followed  immediately  by  its  associated  print  statements 
in  the  program  listing.  This  will  provide  additional  continuity  to  the 
listing.  In  some  cases,  the  print-out  will  follow  later  in  the  program. 
This  provides  output  continuity. 
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3.4 


Sequence  of  Program  Operations 


1)  Reads  number  of  values,  number  of  groups,  data  values,  selected 
critical  chi-square  value,  and  associated  confidence  limit. 

9-13 

2)  Writes  out  number  of  values,  number  of  groups,  and  data  values  as  input. 
19-24 

3)  Data  values  are  sorted  in  ascending  order. 

35-46 

4)  Cumulative  frequency  values  are  calculated  for  each  data  point. 

47-48 

5)  Ranked  data  values  paired  with  cumulative  frequency  values  are  printed. 
49-50 

6)  The  mean  is  calculated. 

57-60 

7)  The  range  is  calculated. 

61 

8)  The  median  is  calculated. 

62-69 

9)  Items  6  to  8  are  printed  out. 

70-71 

10)  Standard  deviation  is  calculated. 

72-75 

11)  Coefficient  of  variation  is  calculated. 

76 
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12)  Items  10  and  11  are  printed  out. 

77-78 

13)  A  histogram  group  width  is  calculated. 

88 

14)  The  number  of  data  points  in  each  group  is  calculated. 

89-109 

15)  The  group  boundaries  are  calculated,  and  along  with  the  group  number 
and  the  number  of  data  points  in  each  group,  is  listed. 

110-117 

16)  Transformation  of  data  values  and  cumulative  frequency  values  to 
functions  of  natural  logarithm. 

127-131 

17)  Determination  of  a  and  3  through  a  linear  curve  fit. 

136-154 

18)  The  mid- value  of  each  group  is  calculated. 

155-157 

19)  The  individual  theoretical  frequency  of  each  data  point  is  calculated. 
161-168 

20)  The  theoretical  frequency  of  each  group  is  calculated. 

169-184 

21)  The  observed  frequency  in  each  group  is  calculated.  Group  numbers, 
observed  and  theoretical  frequencies  for  each  group  are  printed  out. 
185-190 
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22)  Distribution  characteristics  are  printed  out.  These  include  ct,  3, 
data  point  numbers,  data  values,  and  individual  theoretical 
frequencies. 

191-199 

23)  The  value  of  chi-square  is  calculated. 

213-215 

24)  The  confidence  level  for  the  calculated  chi-square  value  is  determined. 
220-255 

25)  The  selected  chi-square,  its  associated  confidence  level  and  the 
calculated  chi-square  with  its  associated  confidence  level  are 
printed  out. 

256-261 

26)  The  calculated  chi-square  value  is  compared  to  the  selected  critical 
value. 

272 

27)  A  statement  concerning  the  relative  sizes  of  the  selected  and 
calculated  chi-square  values  is  printed  out. 

273-281 
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3.5  Program  Variables 

The  following  is  a  summary  of  the  pertinent  variables  in 
the  program,  their  program  symbols,  and  their  output  readings. 


[lajiioJoJit 

?/wQMm 

Output 

V^lgnation 

H&adtng 

no.  of  data  pts. 

N 

No.  of  data  pts. 

no.  of  groups 

N2 

No.  of  groups 

data  values 

array  X 

Data-  as  input 
&  in  ascending  values 

cumulative  frequency 
i/N+1 

array  F 

Paired  with 
ascending  data  values 

Mean 

Q 

mean 

Range 

R 

range 

Median 

DMED 

medi an 

Standard  diviation 

STDDEV 

Std.  Dev. 

Coefficient  of  variation 

CVAR 

CVAR 

Lower  group  limit 

array  B1 

Low  Lim  (Not  Inc) 

Upper  group  limit 

array  B2 

UP  Lim  (Inc) 

Data  pts.  in  group 

array  T 

data  pts.  in  group 

Group  mid- value 

array  VM 

mid- value 

Observed  frequency  in  group 

array  T1 

DBS  Freq  in  group 

calculated  frequency  in  group  array  Fll 

CALC  req  in  group 

Alpha 

A1  pha 

Alpha 

Beta 

Beta 

Beta 

Individual  calculated  freqency  array  P 

CALC  freq 

Selected  critical 

CHCRIT 

Selected  Chi-Square 

Selected  confidence  level 

SIGLEV 

Selected  Conf  Level 

Calculated  Chi-Square 

CHISQ 

CALC  Chi-Square 

Calculated  confidence  level 

Q2 

CALC  conf  level 
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In  addition,  the  following  variables  appear  in  the  program  but  are 


not  output 

array  Z  =  £n(X(I)) 

array  FI  =  -£^(£^(1/(1 -F(X(I)  )) 

=  z  Z  F 
Sp  =  Z  Z 

^  see  part  8  of  program  description 


S3  =  Z  F 


PI,  Q1 ,  and  Q3  represent  the  continued  fraction  expansion 
(see  part  13). 

PROD  represents  the  denominator  in  the  Q(x^/N) 
summation. 

SUM  represents  the  Q(x^/N)  sumnation. 

P2  represents  a  calculated  significance  level  (1-Q2)  which 
can  be  printed  out  of  desired. 
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3.6  Sample  Output 


DLViEixiSlOiM  A(  inu)  rTduf)}  »FiiuO)  .FKluO)  ,p(inu)  »wM(inO) 

_ DiMEi'^s Ji0i4 » pa ( 1  uiD - 

l)I(ViEi>jS10,'j  t  KlUO)  »Hli(xOi)) 

C  nAT«  POINTS  A„l,  the  OlSIRcO  NuMoEH  of  r,„0UPS  hOR  HISTOCRW  INFWWTIOI 


ilPWfillCTrla 


igfSiMgiiBafltg 


c  \/ALuE  are  IiNiPuT. 

/•* 


RtAD(5»?. )  IM 


rEAD(5»1)  U(I) 

.  , r^  /  u  .  /I  \  I  H.-  D  r 


RcAD(5»5)  SIwLEV 


FOKMmT ( I ^ ) 

.V  X  /  T  ^  \ 


FijKMaT  ( Fid  .  0  ) 

T  /  C  1  .'I  \ 


\  A  TE  ( b  ►  i  d  ) 

laf  L.‘  I  T  ^  t  •  ”1  I  i  ivl  _ _ _ - 


w/K  iTE(6M  ?) 

..I  .  .  I  TfcT  (  •  i  K  J  J\lO  _ _ _ 


lAlK  iTE(6»F0i 

.  )  I  TL  f  .  1  1  I 


2  7*  12 


format  (//*  No*  OF  bRolJPS*) 


FORMAT (// ' 


DmTA  -  AS  TiMPUT’) 


C  \/AL0ES.  the  CuMuLmTIVc  FREoJE.viCY  F'OR  EACH  OATa  POINT  IS  CALCULATED. 


iF(x(i)  .01 .  x(i+in  GO  ro  au 


IF  IK  .6T.  1)  GO  10  3U 

TCi  iriii  _ 


RT=X(l+li 

Y  (  I  +  L 


XlI)=RT 

f'/i  Ti"'. 


1)0  800  K  =  lfN 
-Ik  i  .+r 


lin  FjKMrtTl//’  uA  I  A-AsCLNiJIiNiG  ./hLUFj  pAiPr-O  WITH  FKFq,«) 

_ - |iL;-4j^34AlX5A-e-fe^lF:v»4a4i4;-4^iJ^  tj-j-h.  W « "  >  H'  1  ■■> .  b  >  F 1.  U  «-!a-). - - 

b-6^  C 

_ F-  THF . R/W^>£^d£Ja\i.->-4VtEiJ  lMN-t-^UAilDiVri>^-,a£Ul4iX4Jli^JLiUiU-X-0£F-bXCXj:F^^  - 

b;3:.  c  Variation  akE  calluuAI'^u. 

F.,:>.-«£_ - 4; - - - — - - - - - 

5  7i-  iO  =  Ue 

-&F4: _ (4.X_^F — I— T-S-N - 

b  '-3  F  p  0  (g  «  0  +  A  (  1  ) 

- - - - C  i  ■"  O  J  i  4 — - - - —  ■  . .  “  ^  

b  i  *  R  -  A  ( i\! )  “  X  (  1  i 


bx  * 

rj 

■'.  -  •X' 

- 

6  /  •/■ 

— — 
b 

— ly.  t  i.  Je. — 
/  j. 


If  CV  •L.To  1j«}  bO  To 

- - — - - - - - 

0  +  s  « 

_ O^i-X-lbiJ^  A  - - - 

ft . )  To  2 1  i I 

■■Fi-y.|.)— - ( w^i  •)/  fir* - 

i.j  .iLD~X(Lb3 

-.-F-L-U— — ^ — m-4-iTb  iti  »  -  ..  - -  - 

ju f .  i  Tb  (  b  >  2 b 0  3  0  R  f  Oi^iEw 

_ tv4-uO£-V;=-OL*.- _ _ _ _ — . - 

2b5  l4;|»;'i 

c.'-  ^  c  r  i )  n  H  \/  —  ^  I  .  ’J  I  i  F"  V/  ■+■  (  y  {  T  )  *•  >j  ^  ^  V 


S  ruOE  V~Si.gp  I  (  ST UDlV/ /  -T  3 
_ Cx-FiaF.SJ.Puf:  V  /  eg _ _ _ 

WH  i  TE  (  b  »  2'+  u  ) 

_ iuuJJriEl-b-£-2:ib.L...a£D.  JEIv  t  LX'mR _ 

?4n  F'jKMAT  ( // ’  DE»/'. 

-Siu£> _ £4^:V.!.aXjL,^.iA-^-^;.£x5.a.5  i - 

?  2  0  F\;  K  M  A  T  (  /  /  ’  ^V  '  ('J 

^  -1  n _ fcL.jaMi\..I-L5.i<^L>)  El-S-jj-Sul _ _ _ 


CV.aP  •  ) 
RA.jftt 


M'^DTA!'-”  ) 


-T— T.fJF-Q.R.vi ‘.  1  lnM.  CoNLEwMIMb  a  nTx.T.;ftR/L4 .  b Q.,  .iHj^.ijAxA. . Ib^aBLiPijrEil.  _IWi5 _ 

C  le-JCLDOc.b  UlVluTigG  DmTh  x'liO  (3kOi.JP;3  >  :5Fi_FuTif'ie3  1  He  HTU-VALUF  QF  TrlE 


<i 

H 

_ C . fiR  OUK  >.i,iMU.-..UU..UiM-Li.:-Ii.i— 

c 

s  fh  p-x  t  N  )  /  ,u2 

K  •  II-  III'.  |A 

. ! 

1  ,b— J  Pi.  ..b./AL.,a,...A3L: 

B  H* 

i)0  40  0  K4l(.M2 

/ 1  n  i~l  _ 'f  i  K  >  ell  - 

9i=^  J^i 

- .Sj^ir... . .  K— X  - - - — - - - - - - - — 

lOn^ 
lU  /* 


ih,  332  I-1»'\!2 

-J-U-tiEKF^J-iilii-iTcL^T  X-,(-iJ _ 

Ir^fiiTpr^t-  oLlo  NJ  uO  TO  .*534 
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1  1  5* 
-1  1  O  *- 

117* 


12^* 

-12*i^ 

12".* 

- 1 2  n  *— 

12  f* 
-12.H*- 
124* 

131* 

152*- 

133* 

1 


I  )u  52 I  —  1 »  N2 

../Tl— ■ 


_ ^444:EX^*4a4jJ----I-4ilX-(-x^  - - - - 

S2n  C'jNTlf'lUE 

r.tjn  PjKM'^T  L/y.J _ _ - - 

■tz  >,»oup  uo„  LI:,  UP  LiM  PT^M 

Ri it  (-  OKMaT  (  0  X  *  l5 »  2tlb »  I"  12* s) 

r  THF~^A  1 A  vALIItE  Awn  Cu'''L(LMTiVH  FRfOUF  .CY  VALUF';  AKF  TKANEFORMEO 

r  Jf)  H R  '  i S F n  1-tJ A L-l MRAk  i  Iii'l  OF— liji- — L YRE — Lu  ( l-N  (  1  / .(.1  ~F  .(-X.1  — I 

C  LINFAiv  FUwCTIo-J  of  L-MCY)* 

_C - — - - - - - 

R  i  tP=  <  t  N  )  /i''i2 

_ ij^j-  .b-n5.  -I  - 1  *  - - - — - — — — 

Rnb  ^ I  1 ) zALOw ( A (  f ) J 

_ I).  ,  aia  .K,::.!  i .W _  _ _ _ _ _ _ _ — - - — - 

Fin  F  i  I  K  )  -  ALuO  I  AlOu  ( 1  •  /  ( 1  •  —E  i  '<')))  ) 

c  THE  Values  oF  alPmA  AwO  FOR  P  (  X  )  =1-1 /EXp  (  I  V/p,ETa  )  **ALPHa  )  ) 

r  awP  r  ,a:  PI  ii  ATFii  ■AuS*^u-.A--J-Eu*SJL..i>-CYLlAilEb-J--J.-l  « - - - - 


■ 


14.1*  F2n 


14/*  53ll 


144* 


bP-F 1 ( I ) +5o 

■  I.-.  bPQ  _ _ _ _ 

Si4  =  Z  I  I  )  **2  +  44. 

JiUEH  AP-iJliiil.  =E2?..*  SaI  /Aii*a'}-S2.**2±.. 
Hiit-P  (  S4*s3-S2*S1  )  /  i  I  -.2**2  ) 


197* 

199* 

..-tdXlii* — 

9U.-i* 

- - 

iiO  T* 

— gjtiH* - 

/* 

'  j  n  .1 


d:G9* 
-dJM^ — 

^16* 


lJy>  7H'1  I-lP'J 

,Ai  Pi  1  )- JLiXiXijt-tL(J^l 

COimTI ''JuE 

-ti.  iiifidiiXlZJ _ ixSilUlT _ _ 

F  jKMAT  (  ’  {>jO  •  V/  A 


v/AlDE 


I i'j  bP Ol  )P 


Calt  PR  run 
IM  ■''.ROUPM 


fgkmaT { // ’ 
M  A.  T  ( /  * — 
FwKMAT ( ’ 


uT:3Tn(Tbll  I  TuN  Ci-iAkACTh  RiS  i  TCS  »  ) 
Al_F^^A  tiElA*) 


PbP  F  jKMATI/*  HaTm  DhTA  CAi,C») 

C)  <,|  I  .  uiV/i.iT  (  S  D  i  wA.  II^^' 


CAi,C»  ) 


7d2  F'oKMMT(0Xf  Ab»?t.lt).R) 

C  THE  CHi-S9lJAf7ti  VALUE  E  :^K  boOO,MESS-OF -t- 1  I  TEST  CALCULATEU. 

C  >  i  1 S  iJ — 0  0 


HHnMH 

Abl 

r 

_ i  'Xi»t — QisJ _ Lgy  J.  _ 

CniS(>)=CHlSi^+(l./lN+l.  )-P(I)  )**^/H(i) 

21  !* 

C  THE 
r 

Co.^FILENCE  lEwEi.  ASbOClMltn  with  fHt.  oALCiiLATcn  CHT-EoiiAPE 

1 F  1  L  1 1 F  T  F"  tif  M  1  f\l  1-  n  . 

El  9* 
.-Iti'Zt.'i* — 

r 

— E^;-;;c.v.a-b..Q - - - - - 

“ 

221* 

22d* 


Y-SQKT ( YR ) 

I  i  i  -  (  “  1  o  )  *  *  i  \l 
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2.bi*  1  j-»..!UU  =  PKO.»*£' .  »K 

_ l-iJZ _ s>xji'.i=f»LiM±.lt  2»  ^l/PRu3. - 

,^.bd*  ( 1  • +Sl)  J /f- XH  (  I  •  ) 


1  no 

H2=l *-^2 
►<•  l  T  F  (  hi »  /  S  ) 

^  fcx  ^  * 

-i;  i  TE  (  b  »  /  )  » 

1  TP  (  A  »  /■  r  , )  THOk  T  1  •  ST  jLi-  '/ 

mJ  L  >  .-Y. 

db:)^ 

w  -  1 T  F  (  b »  /  b  H  ) 

ut  •  1  T^'  (  K  »  7  ) 

dbl* 

n-k.? 

wR  I TF.  (  6  >  /  Sb  ) 

P.  ik  VluT  (  /  * 

CnlbOf  Qir; 

Eel.  eC  rE:j 

5i-  llC  I  cy  * ) 

dn.'i* 

7d4 

fV/RMAT  (  * 

Cril-S-ollARe 

CuNF  lEvTu’) 

7*^  A 

2hb* 
jiiS  (;>*■■ — 

7rt4 

_ZdG) — 

F>jK  ViaT  (  /  ’ 

_ Fi,.h!MAT,l„* - 

CaLv, 

—  Ehj!  IhI.'L-I - 

Call  » ) 

ruMi-  eFVrLM _ 

78rt  FOKMaT  (  Ox  »  c  Is  •  b »  r  I  •  b ) 

- G - - — - - - — - - - - - 

dh':)*  C  THF  CaLGULATEl)  LHI-SQUAkF  VmL.F  IS  COi..PARc.n  To  I  He  SEuPCTpH  VALUE. 

r  A  c  i  A  I*  iL  iVi  c*  Ki  T  R  i- iA  P  i  j  T  i\i(^  Tut-  C  v)  M  P  A  P  T  s  O  i  M  I  S  H  i  ^  Tj\i  T  rL* _ 


27  /* 

'J~7  4  yHf 

“7  An 

GO  TO  7Su 

u,  fV  1  T*-"  f  I'-s  e  /  A  >’  ) 

2  /  )* 

Fl  t  1  X- 

W  <  ITE  /"bA) 

.,1.4  i  T \r  1  9  7  .'Ar^  ) 

2hi* 

7b0 

.  X-  1  Vj  ■  ■'  ■*  ^  *'  •  ■  ■  '  ■"  ' — “ 

continue 

fil 

- 5^4-U - k.-ih^iaT  1//4 — Xi:A£„rA!,x.iii, ft  I  FL.L._t/.Ai-iir  ri-iT-<^.oi iapf  tc;  t  t-c;Q  tm/ui  opO 

74?  F0KMAT(*  EisiljAL  To  IHF  CRjTICaL  V/aLuE#  I  Hf  HYPOTHESIS  THAT  THE’) 

— - 74-4 - M A T  (  * — ClEa74-UE-0 — r  i, Bt >T X Oig  t--iii->C.l  Ti7?J T c,  FOiin  OaMixIOt  _  _ 

746  FoH;MAT(*  Pc.  JEC  I  Eo  ,  *  ) 

r:*  +- - ?,>? - — T44.E  - .  C-A4-UiTo.a  i  F,  j  t/AUJE — liF. .  C  M  17?_S/.vl  I A  P  t:.  TS  OP^ATFP  THAMf) 

?!]'/*■  764  FuWMmTC*  TnE  CPll'ToAo  v/At_UE.  ThF  HyPOthE^TS  IT-IAT  THE  OERIVtO’) 

— =:)4A44 - 7-64 - EA^-MaT-C.^ — uLi-S-TP  IRUT  i  Di  1  |-l  ii  irTIO^l oOiiD  MijQT  hF  Bp.  li- Ttfo  .  t  ) 

r‘,  ^  f“  I  ;  U 


// 
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.■i7(i0i)+(l5  ,37ni.n  +  Uh 

4  30  01) -fas - >AbXlnH-*-Q.h- 

b0o00  +  n5  .SlOitO  +  Oh 


•  oauOil  +  Ub 
« u  +  l)5 

•  r>4iJOo  +  US 
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CHAPTER  IV 


TENSION  OF  AN  ANISOTROPIC  PLATE  WITH  ELLIPTIC  CUTOUT 

4.1  Introduction 

Consider  an  infinite  two-dimensional  plate  with  an  elliptic 
cutout,  subjected  to  uniform  tension  at  an  angle  (j).  The  dimensions 
of  the  cutout  and  manner  of  loading  are  shown  in  Figure  8.  The  solution 
may  be  obtained  by  superposition  of  the  solution  to  two  problems. 

The  first  problem  is  an  infinite  plate  with  no  cutout  and  tension  at 
an  angle  <!);  the  second  problem  is  a  plate  with  an  elliptical  cutout 
with  no  loading  at  infinity,  but  tractions  applied  along  the  cutout 
surface  corresponding  to  a  uniform  stress  at  infinity.  These  two 
solutions  combine  to  give  a  stress-free  elliptical  boundary  (see  Figure  9). 

4.2  Problem  Solution:  First  Problem 

The  solution  to  the  first  problem  is  a  uniform  stress  given 
by  equations  (36)  below. 

a  =  p  COs2(j) 

X 

a  =  p  sin^t})  (36) 

'  y 

T  =  p  sincj)  cos(|) 

xy 

pKoblzm  Solailon:  Second  PAoblan 

The  second  problem  will  be  solved  using  the  stress  function 
method  of  Lekhnitskii  [15].  Stresses  and  displacements  are  obtained 
in  terms  of  the  stress  function  '^’|^(Z|^)  and  its  derivatives  with 
respect  to  Zj^. 
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(37) 


X 

=  2Re[${  +  $2] 

=  -2Re[yi${  +  \i2^2^ 

U  =  2R^[pi$l  +  P2'^’2] 

V  =  2Re.[qi$i  +  q2®2] 

The  parameters  p  and  q  of  equations  (38)  are  functions  only 
of  the  material  compliances  g.  .  and  the  roots  of  the  characteristic 
equation  for  the  material. 

-  26x61^^  +  (23x2  +  366)^^  -  2326U  +  322  =  0  (39) 


Pi  = 

3iiPi^ 

+ 

3i2  ■  3x6V^1 

P2  “ 

3iiP2^ 

+ 

3x2  “  3x6>^2 

Pi  = 

312^1 

+ 

322/1^1  “  326 

92  = 

3x2^2 

+ 

^Zllvz  ~  326 

(40) 


The  constitutive  relations  employed  assume  only  that  the  materi¬ 
al  is  midplane  symmetric;  otherwise  the  material  may  be  generally 
anisotropic.  The  roots  of  the  characteristic  equation  are  in  general 
obtained  numerically,  as,  for  example  in  reference  [16]. 

The  Z|^  coordinate  frame  results  from  a  linear  transformation 
from  x-y  space,  where  px  and  P2  are  solutions  to  the  characteristic 
equation. 


Zi  = 
Z2  = 


X  +  PxY 
X  +  P2y 


(41) 


The  stress  function  may  be  represented  as  an  infinite  series 
in  the  variable  as: 


oo 

\  logc  ^  +  E  5,^ 

The  variable  is  the  result  of  a  conformal  mapping  of 
such  that  the  elliptic  boundary  maps  onto  the  unit  circle.  The  mapping 
is  given  by: 

*  ‘^k  ■ ■  ’■*'k‘'> 


The  problem  is  now  reduced  to  the  determination  of  the  stress 
function  Since  the  net  load  on  the  interior  boundary  is  zero  the 
coefficients  Aj^  reduce  to  zero.  The  remaining  coefficients  are 
determined  from  the  boundary  tractions  on  the  cutout,  which  may  be 
represented  by  a  Laurent  series. 

The  boundary  stresses  are  related  to  the  surface  tractions  by 

T  =  a*ri  +  T  *n 
X  x"x  xy  y 

T  =  T  *n  +  a*n 
y  xy''x  '^y  y 

The  components  of  the  unit  normal  to  the  cutout  surface  must 
now  be  determined.  The  equation  of  the  ellipse  may  be  given  in 
parametric  form  as: 

X  =  a  cose 

y  =  b  si  ne 

Consider  an  incremental  arc  length  dS  along  the  elliptic 


surface  (see  Figure  11). 


■y 

Figure  11:  Local  Coordinates 


68 


dx  =  -a  sine  do 
dy  =  b  COSO  do 


=  dy/dS 

n  =  -dx/dS 
y 

Tho  strGSSGS  a  *,  a  *  and  at  thG  boundary  of  thG  cutout 
X  y  xy  S’ 

arG  known,  sincG  thoy  must  corrospond  to  a  uniform  stross  appliGd  at 

infinity,  so  from  Gquation  (36); 

a  *  =  -p  coscf) 

X 

=  -p  sincj) 


*  =  _ 


p  sin(f)  coS({) 


ThG  boundary  tractions  arc  now  fully  dotorminod: 


^x  X  dS  xy  dS 


T  =  _  *  ^  _ 

'y  ^xy  dS 


y  ds 


ThG  boundary  tractions  can  bG  shov/n  to  have  thG  following 
relationship  with  the  coefficients  of  the  stress  functions. 

/ (T,  t  1Ty)  dS  =  -(P,  +  iPy)  !;■ 


+  L  [(yi-i)  +  (v2-i)  Ap  ] 


2m-^  e 


+  Z  [(yi-i)  A  +  (y2-i)  A2  ] 
n=l 
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Representing  the  integral  above  in  terms  of  a  Laurent  series, 
it  is  possible  to  solve  explicitly  for  the  coefficients  needed  for 
the  stress  function. 


/(Tx  t  iT^jdS  = 


B  e'"®* 
n 


E  D 

m=l 


■im0 


m 


(51) 


1+i^^ 

.2{yi-y2) 


+ 


I4di!2_  D 

2(yi-M2)  ni 


1+i^^  n  +  Q 

2(y2-y  1 )  m  2(y2-yi)  m 


(52) 


Returning  to  the  expressions  for  the  boundary  tractions  (49) 
the  value  of  the  integral  in  (51)  is  computed. 

*  ^ 


/(Tx+iTyjcIS  =  “'Vf) 

=  f  fix  *  +  ia  *)  a  sine  d0  +  (a„*  +  ix^,,*)  b  cos0  do 
J  x-y  y  X  xy  ' 

=  -i^xy*  ^  ^  ''^xy*^  ^ 


r  I 

y(T^+iTy)dS  =  1^ 


[■ 


T  *(a+b) 

xy 


+  1 


(aoy*  +  bg^-^) 


(aa  *  “  bOxx*) 
'■  y  X  ' 


10 


-1  0 


(53) 


The  coefficients  of  the  Laurent  series  of  (51)  are  now  easily 
calculated,  and  lead  directly  to  the  coefficients  of  the  stress  function. 
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“  (“^ii^i  )  / 


$2  ■■  {~^2l^2  )  / 


’  The  equations  (37)  and  (38)  may  now  be  used  to  obtain  stresses  and 
displacements  at  any  point  in  the  plate,  due  to  a  given  load  p. 

A  computer  program  has  been  written  to  implement  this  solution, 
and  a  listing  appears  in  the  next  section. 
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4.3  Program  Listing 


C  stress 'FUKicTlOM  SOLUTION  ^  ElLlpTjC  CUTOU?  !N  jNr|NlTE  PLaTe  -  jaN  7;  ' 

C  AA  s  semi-major..  Alls  _  I  X’^AxiS  ) 

C  r  B  «  SEMI  *  minor  axis  s  Y  axis  ) 

C  ThETao  ■  lUfLlMAT  lON  ANGLEiriEG)  V.  jTh  X»AXiS  OF  TeNSIOm  a?  INFIUJTY 
C  PP  ■  VALUE  OF  TENSION  Aj  INflNlTV 
C  B£TA,.S . FLEilBILlTy . ,mATR  |  x  Fnf'  M  ATrR  t  AL 

C  MU  E  SOLUTIONS  TO  CH A R A c T E R | ST ! r  fOUaTION  FOR  MATERIAL  '  r 

COMMON  /  HAP  /  ,.Pi2  5  sQL2).»BETa(6J  sMUC|)  ,  .m 

COMMON  /  GEOM  /  AAj  B89  xCOR  f  5s  J  0  VCOR  « S5  3  <.1  i, 

...  DlME,MSJ.mx__TlTLFL.2Q)  iX(2Q0!  tY(2Q0HGQ|S&S 

complex  mU,ZFTa(2I sZ«2? jtSlGXtCSIGYsCTAUsCU^CVsRADCL/ri sOPhUJ) 9  ,  f 

,l...PHai.2,UA  (2  9a-)  *J,.»P*PaO . 

REAP  (5,303)  Nc 

14  .N £..._=  t  . . . 


IF(NCiLT»0)  STOP 

Nc . -  . - . 

RE  AD  (5,  ion  title 

.M|TXiA.».iQ,l  )  .T-ITLE . 

reao{S,io2) (bCta( n  ,  s-i^a) 

M4t.e.ia-,-2DU  (BrTA(  n  ,  t=>i  ssi 
REaD(S,jP3)  (Mint)  ,  IB?  5,25 
WRiTEi4..202)  (Mud)  t  1  =  1  «2) 

READ(5,to3>  AA,  B8  ,  PP,  tNETaO 

.  w.,R  1  T£li. ,1  2.Q  2 i.A,.A  , . fi.B ... .  P P  ,  t  H  e  T  a  0 

theta  =  3t 141592653469  •  ThEtAh  /  180*0 

.  IJ  ..  S-lQ-t-di  I  i  0  ) . - . 

calculate  uniform  stress  FfELO 

S16  1..  .■.,.9,_PP.«.  C..05(  theta  )  «»  CQs(tHETA) 

SI62  *  •  PP»  SiN(THETA)  e  SIm(THETA) 

S  .1  G  j  2  ,.P  p  ? . S  J.  N  .(THETA)  «  CoSiTHETA) 

W  R  I  T  E  {  6,  2  9  9  )  S  1  G  I  1  S  I  G  2  ,  S  1  G  1  2 
299  F.OhMaT(  lOX.3E20»iO  ) 
calculate  A(K,N) 

B....  .  SIG!  2  «JBB“AA»  Tj  ♦  «5  ».(aA«S 

P  B-.S  »  SIG12  •  (BB'S'AA)  »  Tt  *  e5  ®  (aA»S 

„  A  (.1,1  ) .,  « J  J  *  t,  -T  I  •  MU  (  2  )  )  /  <  2  *  ®  ^HU  (  3  * 

1  <*.  (I,  w  T1  •  MU(2)  )  /  (?,  *  (Mud!  -  MU 

A(2ilA  a  (1.  <■  T!  «  MU(li  )  /  (2*  «  fMu(2) 
I  ♦  (I*  -  TS  *  MUCd  )  /  (2.  «  {MUC2)  "  Mu( 
C AlCUL ATE  P tJ  J  ANp  0(1) 
calico M ST 

generate  .PES  IREO  interior  SqLUT  {On  PO  !  NTS 
CALL  GEM(QQ) 

D  O . 14  ,  J„.B_„l  tS5 

Y ( J  >  “  OsO 

3  4  _  X IJ  L .  j«  _  A  A.  *  Q  fi  t  J 1 . . . 

34  X(J)  »  qQ(J)  ♦  0»227 

. .  .. 

J  «  JJ  T  SS 
X  (  J  )  ■  n  f  0  .  _ _ _ _ 

35  Y(J>  =  Qg(jJ)  ^  1»02S 
3  5  ,  Y  (  J  )  . 8  0  .  ■*■  0  Q  (  J.J  I. 

write (6, 204  ) 

DO  39  J  9  1 1 dn 
DO  *4  0  K  B  1,2 
KEY  «  0 

SiGf'falgO  79 


»,(AA«SIS2  +8B*SI&1) 


Ml  ■  .7 
''0 

£t.PK-‘X'U 
El.P  5  UO'/O 
g;  L  F  s  0  g  *?  S 
r  5...  p  ;i  u  3  0  0 
F  L  7  I  0  I  G  5 
E  L  P  3  0  U  Q 
E  L  P  5  TM  4  5 

O, i.  '4  j  2  5 
ILPI  0  3  3Q 
ELFIOiSG 
liLPIGl  ‘'U 


E  I-  (  G  3  (  G 
E  L  P  I  '  'J  n 
EG  PI  ■'  i-u 


B  B  «  s  1  r,  1  ) 

£  L  P  4  0  i  /  0 

;:U  (  2  1  )  )  «  B 

F  L  P  1  0  i  7  B 

)  )  •  D 

£  L  P  1  U  i  8  U 

M  U  (  1  )  )  )  *  B 

.£ k p  |C3 1  a 5 

)  •  D 

EkPlOlfQ 

£lF id  j  9S 
E  L  P  1  0  2  U  u 
ELP I02Q5 
ELP 1 02 1 Q 
EkPIOSii 
ELPI3  k'G 

E  L  P  I  G  '  i" 
g;  L  P  I  0  M'l . , 
■€tP|  0  M'G 

elf  10?  > 

.EiPinp  M 
LLP  I 0  2SG 
£tP  I02SS 
g  t.  P  J  0  2  6  C 
ELPI02GG 
£  L  P  I  0  2  7  G 
E  L  P  I  0  7  /  s 
E  L  P  i  0  2  8  U 
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-  — . . .  .  . . . . .  .  . -  elP»02B5 

RADCL  ■  Z(K)  •  ZJK)  *.  AA  •  Aa  -  MU(K4;  •  MU(k)  •  BB  *  RB  EUP|0290 

2ET^riC^-*-4  Z(i4— *  sign  ♦  cSaRTtKAOCU  )  /|^A  •  TI  »  *  frS^  )  EUPIO29S 


IF  (key  *E9i  I  )  (50  TO  50 
RZET4  «-f?EAi(Zf4MKI 
TIZETA  •  Alf^AGf  ZETA(K)  ) 

. T M  --.- Zf^-A-JA-A-Z-tT-A- A-  T  I-2-f  T  A--*-  T  I  2  E  T  A  —  . . . . .  . .  '  " 

ZEE  *  Si3RT(TM) 

. . --  ■  .  . . .  - .  - . .  . 

IF  (ZM  ,UT*  ItE'S)  GO  TO  SO 

SIGN— ■-»-Si6N-  . - .  -  . . -  . . 

KEY  ■  I 

. ...... . . . . . . . . . - . — . . .  — . . 

50  continue  y 

PJilJja__*__AlKUl>^ZfTAtKl-  -  . . . -  . . -  — 

MO  OPhI(K)  •  *•  A(K,l)  /  (SIGN  •  CsQRT(RADcL)  •  ZeTA(K)  ) 

CALCUL.41£_ixR£SS£$^..AN0  JJISPlACE-mENTS . 

CSIGX  «  MUIJ)  *  MU(l)  •  OPHl(l)  ♦  MUI2)  •  Mu(2J  •  0PhI(2) 

CS4-«-Y^-»-oBH-U-4  >■ . .*.-.4JP>i.M.2 1 - - -  - -  . . . . .  . 

CTau  ■  MU(t1  •  DPHKI)  ♦  Mu(2>  •  DPHI(2) 

XUJ~  jlP444  *-P^f444-^ ♦  F(  2I  -  - 

CV  8  Qd)  •  PHid)  ♦  0(2)  •  PHH2) 

•  PEAt  t  CSlGX  y-/Pp  »  SIGl 
S16Y  =  2«  *  REAU(  CSIGY  )  /  PP  »  SI62 

.  ,  . -.*._S16.-l2  . - .  •  — . -  - - - 

U  ■  2.  •  REAK  CU  > 

V  .a  Z  .  -A  REAi  t  -CA(-  ) 

WRITE(6,20S)  X(J)»Y(j),SlGX»SlGY»TAUtU,V 

39  COfil-I-MUE _  — . . . -  . . . .  - 

60  TO  IM 

101,.  FORM, . . — . . . .  ■  . . . 

102  format (AE|3t7  ) 

1 03  F.ORHATj.  MEZOf  iO  )  . -,  -  . - . . 

201  FcirmaT(  1X,6E17.7J 

202  FORMAT.(_^l  A,i.M£2axI0  >  .  . - .  . .  . 

20M  F0RMAT(///////  6Xi*Xi,9Xi*Y*.10X,’SIGX*,9Xi*SigY»|9X,*TAU»»11Xi 


EUP10300 
EI.PI0305 
Ei.P|03l0 
ELP10315 
ELPI0320 
EI.P10325 
tUPl0330 
E1.P1033S 
EI.P103M0 
EI,P103^5  - 
ELP10350 
EI.PIO355 
E!,P10360 
EUP10365 
ELP10370 
tl*Pi037S 
EI.P10380 
£LP  10365 
ELPI0390 
eLPl0395 
EUPIOMOO 
tLPlOHOS 
Et.PJ0M10 
E(.P10415 
ELP10*i20 
-EUP-10M2S 
EUP10M3Q 
EI.P10M35 

elpiomho 

eepiomhs 

EUPIOMSO 

elpiqmss 

EUP10M60 


.  1 . .  » UJLaZXa,.».y.L..  JlA-J.,  . . . . . . -  . . . . . . . . . .  .  .£1,6^104.6,^ 

20s  F0RMAT(2Fl0*9f3Fl3*2|2FlM.8  )  EUP10M70 

301  EDRMATI  .J3JL _  -  - . -  eUPiOHJS 

302  F0«MAT( iH|  ,13)  ELP10M80 

EN|l_  __  „  . . - . . .  .  El-P  10485 
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SUrPOUT  I  nE”  gen  (  Q  )  ' . 

cdordimAtes  for  euliptiC  Cutout  solution  pojnts 

common  /  GEOM  /  AA,  B8,  X CoP f B5 ?  s YCOR ( 55  ) 

d!H|:NSI0M___QJ.55l 

CWfiG  ■  .COl 

Q  U  )  .  ■ D  •  0 . 

DO  50  J  «  2|SB 
IF  (U«!6T,  I  II  CHmG  «  f0I 
ir(J.6T,20l  CHnG  »  .10 
IF ( J,fiT,29)  CHnG  *  *50 

|F(J.GT*<46)  CHfiG  »  1.0 

]F  U...S.Lt5G.) . ChNG,.»  5.0  , 

N  a  U  •  i 


elf  I  5000 
ELP150Q5 
E  i,  P  I  5  Q  1  0 
L  L  P  I  E  Q  I  S 
E  L  P  I  5  0  2  0 
£LFI5D2i  ’ 
£  L  P  1  5  Q  i  0 
ELPiSoSS 
ELP ISOPO 
gLPiSD'ii. 
ELP 15050 
eLPl5D5&,. 
OLPibO&D 


5  0  Q(J)  »---C!iN.)..A.xHN_G..  -  ELP|5Di)5 

return  ELP;-'1^0 

END  .  -ELrUVS 


subroutine  const  EUP20000 

COf-MON..._/.-.HAD.  /.-,,E.U)  .Qt2)  .BeT-AU*  .MU«2)  FLP;;000& 

CD’NPLEX  P,Q,MU  |LP2Q010 

DO  4  -K  •..-}-|2  .  .  .  .  ELf'200|5 

P(K>  ■  BfTA(I)  *  MU(K)  •  MLKkJ  ♦  BETAJ2I  -  BE?A(3)  *  MU(K)  EUP20Q2Q 

4  Qfl^5 )  • MULK  )  -■  BET  a  MI  /  Hy  C  K  )  -  BE  T  A  (  5  )  .  ELP2QQ2S 

ELPFOn^O 

. .  . . . - . — £UP2003.i 
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